IDENTIFICATION OF SHIP MODEL 
MOTION PARAMETERS 



by 



Jorge Manuel Delgado Beirao Reis 



IDENTIFICATION OF SHIP POD EL 



MOTION PARAMETERS 

ty 

JORGE MANUEL DELGADO BEIRAO REIS 
B. S., ESC OLA NAVAL (PORTUGAL) 

(1966) 

SUBMITTED IN PARTIAL FULFILLMENT 
OF THE REQUIREMENTS FOR THE 
DEGREE OF NAVAL ENGINEER 
AND THE DEGREE OF MASTER 
OF SCIENCE IN 
MECHANICAL ENGINEERING 
at the 

MASSACHUSETTS INSTITUTE OF 
TECHNOLOGY 



June 19?1 



LIBRARY 

MAVAL POSTGRA£<IATS SCHOOB 
MONTEREY, CALIF. 93.94Q • ^2~ 



IDENTIFICATION OP SHIP MODEL MOTION PARAMETERS 

By 

JORGE MANUEL DELGADO BEIRAO REIS 



Sudinittsd to the Department of Naval Architecture and 
Marine Engineering and the Department of Mechanical 
Engineering in partial fulfillment of the requirem:ents for 
the degrees of Naval Engineer and Master of Science in 
Mechanical Engineering 



ABSTRACT 



The applicability of the methods of systera identlfl” 
cation to the determination of ship model motion parameters 
is studied taking as an example a miOdel of the Mai’iner 
class hull form. 

A mathematical expression desex^ibing the dynamical 
behaviour of an ocean vehicle is first presented. Then the 
values of the ship model motion parameters to be identi- 
fied are predicted using available techniques. Data 
representing the dynamical behaviour of the hull form 
under controlled conditions is obtained v/ith model tests 
in the tov/ing tank. This data is tPien used togethex'' with 
the fii’st estimated values of the motion partxrae tens to 
determine the corrected values of tPioae pax^ame tens . 

The sensitivity of convex'gence to sevei’al noise 
factors caracteiistic of the process is also studied. 

The results indicate a faii’ly sti^ong dependence on 
those noise factors and on the initial uncertainty that 
in general decide the final convex'gence . This fact seems 
to indicate that the particular technique used and the 
particular input adopted are not in conjunction the most 
appx’opriato for this type of dynamical behavioux"*. 
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CHAPTER 1 



INTRODUCTION 

Consider a system whose inputs are known and the output 
ohservable. System ident if icat ion , in general, has the 
objective of determining the structure of that system 
using the known inputs and measuring the outputs in order 
that one may be able to predict the beha.viour of the 
system (outputs) to the va.rious excitations expected from 
the environment (inputs). One such a system of interest 
in naval architecture is the ocean vehicle. The naval 
architect wants to be able to predict the velocities, 
accelerations and other important motion responses that a 
ship wi.ll experiment at sea in order to effectively design 
■^or proper response, structural integrity, and other 
critical aspects of operability. 

The id en Ixif ics t ion of a system of completely unknown 
structu.re is, in most cases, impossible, least for the 
tine being. There are, however, some types of systems for 
v.'hioh identification is possible. Those consist of para- 
meter identification problems where the genera] structure 
of system is known but the particular values for a set of 
param.eters are to be determined. The ocaan vehicle ’■^as , 
fortunately, those characteristics. 'Using rigid hody 
dynam.ics, one is able to develop t'^e so-called equa.tions 
of motion which descrihe the dynamical hehaviour of the 
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vehicle. Tn these equations there are some parameters 
(ship motion parameterr or hyhro'i;^ma’-iic coefficients) which 
are not known or at host v/hose value is known with a certain 
degree of uncertainty. It will be necessary then, in order 
to fully predict the behaviour of ocean vehicles at sea., 
to obtain the correct \’’alues of those parameters or, at 
least, by selected nea.sursnents , to d '^crease the degree of 
ij.ncertainty . Unf ortunately these measurements are, more 
often than not, combined v.'ith noise from th“ instrumients 
themselves or from the environment at the tine the measure- 
Ficnts take place or from both. These types of noise need 
to be t'^kon into accm^nt when identifying the parameters. 

An example of the fo'^^'^going is the case of experiments in 
a towing tank to obtain data from models to apply to full 
scale vehicles. Measurement noise exists whenever one 
uses any type of transducers that always show some 
hysteresis and are characterized by .a specific dynamical 
behaviour. Environment noise will be cav,sed among others 
by interference between the model .and undesired v.’aves , 
and vih^a^tion of test er[uipm.cnt, Of course these e.ffects 
could h""’ included in the ma.thena tical m^odel of the vehicle 
but it would make the analysis of its behaviour much m.ore 
dif.ficult. Thus it is seen that sometimes it’s preferable 
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to a sinpler model ■Tor the 'behaviour of tie system 

and introduce the effects not accounted for in the mathe- 
matical model as noise (environmental noise) whose general 
nature is assumed known. 

Many of the available techniques of system identifi- 
cation were developed for a dynamical system whose 
behaviour may be expressed in a so-called state space 
representationfl]] ; 

X = f(x»u,t) ; 2l(tQ) = Xq (1.1) 

11 = ( 1 .?-) 
where the vector x expresses the state of the system at 
ti"'e t ^t;> * "^he vector represents I'he cu.tput o^ tlie 

system at tim.e t and the vector u. the inrut at ti’i^e t, 

Xq is the state of th'^ system at time t^ (initial conditions) 
and X is the derivative of x v/ith respect to time'. f and 
hare vectors that describe the structure of the system. 

The main assumption to be made at this stags about the 
above system in order that meaningful results ce.n be 
obta.inod is that there will be unique values of x(l) 
v(t) for a given u(''') and x(lo) “ — o ^ '.hat 

t>to. 

The purpose of any technique of system id ent i "'ioat ion 
'Will then be, by operating on u(t) and y(t), the deter- 
mination of f e.nd h. As said before it wall b-'' , in 
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general, not possible to obtain the eolution to this pro- 
blem, Therefore, the 9.seumption of a. known strooture 
except for a set of unknov/n parameters is made and the 
system is expressed by the following: 



where the vector p represents the set of parameters to 
be determined , 

In order that the identification of those parameters 
be of some value, it is further assumed that p is time 
invariant in the time interval of interest t^^ t^< t , or 
at lea.'-'t. that its time dcoendence is known meaning that 
once some p(^) , t^<”C^<t^ , is determined, p(t) will 

be known in the interval (]t^, t;^l« 

It must be emphasised at this stage that the tore- 
going does not assume the uniqueness of that representa- 
tion. Only the \miquencss of the input/ou tput relation- 
ship is assumed. That is, there v/ill be 'for a given 
p(t^,t) and v(t), a unique b such that 



It v/ill be possible to sliow the existence of a repre- 
sentation x' , , h' and p’ such that 



X = f (x ,u ,p, t) 
Z = g(x,u,p,t) 



Z(-to) = 




y(t) = t(li(tQ,t) , t) 



( 1 . 5 ) 



y = h'(x'u,p',t) 
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The problem is then redefined as the determination 
of a set of parameters of some state representation of 
the system defined by y(t) = b(u(tQ,t), t). 

Coming back to the case of an ocean vehicle, it will 
be sho’.vn that the resulting equations of motion may bo 
r^^gardod as a state representation of the vehicle 
dyna.mical behaviour and that the assumptions made so far 
are met by that mathem.atical model. 

The foregoing assumes, apart from the unknown set 
of parameters js, no other types of uncertainty. But, as 
said before, it is som.etimes advantageous to select a. 
simn3.er m.athematical model for a given system and Trtro- 
duce the so-called system or environmental noise since- it 
results from interactions between the system and the 
environment. Also consideration must be given to 
measurement noise. A model that takes into account those 
factors will be of the following form 1^2 1 t 



21 I(21»li>]2»ll»‘t) 

'JL " ) 




( 1 . 8 ) 



( 1 . 9 ) 



whore the vector w reoresents the environm-ental noise of 



the system' and the vector y the measurement noise 



"^Errors in the structure of the system will come through 



as environmental noise. 
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Tn addition (case of time varying systems) the parameter 
vector p may have the foilov.’ing structure: 

p = I2('to) = Ro U.IO) 

where £t is known and w’ is noise. 

In order to simplify the above formulation the 
follov/ing assumptions are generally made; 

1. The mathematical m.odel is time ' invariant . 

2. System structural uncertainty will appear as 
an apparent additive noise process, 

3. f'easu.renent noise may also be considered an 
a.dditive noise process. 

, Output mea.surements are lin'^a.r '^uncticns of 
the state ^^f the system and are s tructur?,lly independent 
of the input. 

5. Parameters are coefficients in the model 
structure that are to be identified. 

The system dynamical behaviour may be expressed 
as follows: 

A “ £'2i»y:.»2.) (1.11) 

p - 0 O .12) 

y = H X + V (1.13) 

where H is the measurement matrix, often a dia.gonal 



matrix . 
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The purpose of the ihertif ication v.'ill be the deter- 
mination of o. 

In Chapter 2 the actuations of motion for an ocean 
vehicle will be developed and it v/ill be shown that indeed 
they constitute a state space representation of the ocean 
vehicle's dynamical behaviour. It will be shov/n also that 
they satisfy the foregoing assumptions. 

In Chapter 3 "the methods of predicting the ship 
model motion parameters v/ill be outlined and values of 
those parameters for the mariner class hull form v/ill be 
presented (Appendices 1, 2, 3 )* 

In Chapter the excitation from the environment 
v;ill be described and its several components isolated. 

In Chapter 5 'the experimental procedure followed 
in order to obtain the data will be described. 

In Chapter 6 the different methods of system 
identification v/ill be discussed and those which seem 
more convenient v/ill be used to identify the motion 
parameters of the m^ariner class hull form. 

In Chapter 7 the resvilts obtained in Chapter 6 v/ill 
be discussed and in Chapter 8 the consequent conclusions 
aiid recomm.endat ions v/ill be pointed out. 
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CHAPTER R 



EQUATION OF POTION 



In order to develop the equations of motion for an 



ocean vehicle, it is convenient first to define two 
rectangiilar right handed coordinate s^'ctems, one fixed 
in the surface of undistur'bed water v;ith the positive 



^3 " ^o pointing downward and the second '-'oving 

with the vehicle, the origin of such coordinate system 



being arbitrary (but the ooordinatee (:■' ,y^,e„) 
center of gra.vity of the vehicle in this system 
known) pointing the positive x-^ = x axis to the 



the v'eh, icle ?-nd the x- 



z no in ting do’.vnward i' 



of the 
bo j ng 
bo’.v of 



the 



equilibrium position (see Figure 2,1) 




Figure 2,1 



I 
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The or.iontation of the moving nxis sycten v/ith respec 
to the fixed axis system will expressed by three angles 
(the "modified Euler angles") defined as follov/s, 
assuming that initially both axis systems coincide C3^** 

^ - angle of rotation of the moving axis cystoTf] 
about the z = z^ axis (yaw angle) 

9 - rotation of the moving coordinate system already 
yawed about its y axis (pitch angle) 

<P - rotation of the moving coordinate system 
(already yav;ed and pitched) about its x axis (■>r'oll angle). 

Defining the rotation about the body axis as 
respec ti ve''.y 5/L(about the body X]^ axis), S 5 (‘'^bout the 
body Xp axis) and S^(about the body axis) the relation 
between those rotations and the Euler angles will bo C 3 I: 

~ <^cos9cos'f' - lsin4^ 

Sj, = c^cos9sin'^-r Icos'f' (2,1) 

5^ = -(j)3in9 + 

The motion of a rigid body in six degrees of freedom 
with respect to an avis system fixed in space- when the 
origin is at the center of gra.vity a.nd the moving axes 
are the principal axes of inertia, may be represented by 



Nev.'ton ' s 


ecu a t i ons , na.n e 1 y : 


F 


rt d 


( 






dt 












(2.2) 


11 


= d_ 


( 






dt 




lx 


rn 


i s <a 


ri 1 


agonal na'^rix v;hose elonents are egnal to 
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the mass of tho "body, is tho transl.ationsl A''e^ocity 
”ector of the center of {gravity, of components alon;'^ the 



moving axes , Vq and 'Mq, _F_ is the total vector force 
acting on the body of components along the axes X, Y and 
Z, X diagonal matrix whose elements arc the inertia, 



is the rotational velocity vector about those axes 

• T 



of comnonents o„ , o„ and r_ and M is the total moment 

‘ ij 'o G ~ 

acting on the body of components K, M, and N, Neglecting 
the effects due to rotation of the ear+h in ocean 
vehicle dynamics Abkov.'itz and Tufts [ 5 ] develop 

ecuations (Z. 2 ) respectively with fixed and varying mass 
ar.a nos Tt ion cf center of gravity ^ ov the case of an 
arbitrarv oricin, the first using vector anal y^i^ ard the 
second Lagrange's equations, Therefore only the resulting 
equations will be presented here and for a body of fixed 
3,rid inGr’t'iP. . 



X = mru-fqw-rv-x^_(q'"tr'‘ )+yj^(pq-r) + z^(pr+q)l 
Y = m'”v+ru--pw--yr(’^ +p ’ ) + z.„ (qr-p)-i-x (qp+f ) ] 

^ 'j (j fjr 

Z = r.rwtov-ou“Z_(n^+q^) + x^(rn-q)+y (rq+p) 1 
^ ^ ' G ‘ G * G 



K 



M 



2 



= I p +I„„(q'-nr) + I (ftpq) t I, (q '-r ) + 



(I ~I (v-p V/-!- r ti ) - 1 ’ fi". V ^ ( v: - o\x -f p v ) 

^ zz yy ' G G 

I,,-/ + i,«a+ar) + I (g-ph -i- 

)rm-mx^^(w-q'H-pv) + mz^^(u--rv+qv;) 



zz 



(2.3) 



(2.4) 



M 




1 




# 
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= I (q+rp) + I^^^(p^-q^) + 

• • 

( *f --n( v-pv;+ru) 



(2.4) 



yy XX' ' ‘ c- ' " ' "G 

v.'hfsre unsubscripted velocities and accelerations refer 
to an arbitrary origin and x_ , and z^ are the 
coordinates of the center of gravity. 

Kote that F and M v/ill in general be functions of 
position, velocity and acceleration of the vehicle as v/ell 
of properties of fluid, geometry Gif' the body end of the 
inpu.ts. For a body of given geometry and mass distribution 
and in a given fluid then 






= 




II 


~ 


\ 


[u , V , v^l 


n - 


( — 1 


2. “ 


b-o.J'o'^ol 


f', — 


C + .e.+l 



(2.5) 



Since the purpose of this investigo.tion '.vas the 
study of the applicability of system identification 
techniques to ocean vehicle dyna,mic3 and^ not the 
identification of parameters of a particular vehicle 
(carameters i.mbedded in F and M) it '.va,s decided ai: this 
stage to simnlify the formulation by restraining tb.e 
model to one or n.t nest tv;o degrc'^s of freedom. It 'vill 



!?• 



be only ccnsider^^ henve ( tror’-'l atirn along; the z axis) 
and roll (rotation ahou’’' the x oxie) v/ith zero forv.'ard 
speed. To simplify the formulation even more y» and z 

G 

v/ere '^ade z^^ro. Therefore, the equations of 'motion 

for this restrained mode ( y„-z_=u=u-v-^v^o^c-f=rr=x -v - r=^r.Q) 

L- lx - n ^ n ' 



Z = mw 



(?. 6 ) 




and Z = Z['r,v;, z^ , n,p, , inrut I 
K ^ K[\v,’,v, z^.P.p, <|) , i^-rut I 



(?.?) 



Looking at equations {2,6 and (1.7) one 



o 



ef^s that 





p = d 1 ~ ^ 4 eos9 cos 4* a in 'f' I from Eo . (1.1) 

dt '' dt “ 




see Pig. ?.2 







4 



/ 



Figure 2.2 
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ons 



(^.8) 



Since 'j', 9 and v are zero, the above rolati 
reduce to 

« 

(f) = p 

Zq = w cos (j) 

Ignoring for the noment the presence of w and p as 
arguments of Z and K the complete set of equations in 
state space form would be: 

w =1 Z(w,w, z^,p,p,4», input ) 



rn 



z = w corcb 

o * 



p = 1 K(w,w,z ,p,r,</>, input ) 
I 

XX 



(Z.9) 



(p ~ p 

In order to be able to proceed and apply the existing 
identification techniques to ocean vehicles (specifically 
surface ships), it is necessary to isolate w and -that 
the dependence of Z and K on w and p is linear. Therefore, 
the next step wilJ. be to assume the existence of a 
Taylor's expansion about an equilibrium point, to be 
taken a.s v.’-v:~ z(j>= 0 of the form: 



Z ^ Z(o) + oZ_ 
6'.V 



6 inn 



+ 



5Z 

5'"/ 



V' 



6Z 

67 



5Z 

'O -I- 6 .;s 



fi -5 !- 



input .1 6^2 


+ 1 if z 




2 ! 6vy' 


0 





(-^. 10 ) 



0 
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1 


6^Z 


• 1 

4* - 


6^Z 


t 

v; z + 
0 

0 


1 


6^Z_ 




9 ! 


6W6w 


CV 

o 


6w6 z 

r. 


2 ! 


6w6 n 


nd 


.in 


the same v.'ay 


f or K 


• 







( 2 ,] 0 ) 



It is Tcriov/n fron potential theory f 5 1 that for a 



body moving in an infinite fluid av-^a^^ from the surface 



there are no dependence of force or moments on ]:)roducts 
of accelerations with velocities or accelerations. 

It will be necessary to accept this assumption for 
surface vehicles without any other justification, and 



neglect also all other terms where acceleration is 



multiplying any other variable. 

With the vehicle ?n cnui librium also Z(o) = 0, 
further simrlif ice tion can. be achieved if v-’e 
consider eo.ch elemjent of the expansion, which will be 
considered in the next chapter. It will be enough to 
say for the moment that due to port and starboard symmetry 
(both geometrically and inertia 1) the vertical motion 



of the vehicle does not produce any rolling moment, at 
least to first order, and therefore those term.s ca.r be 
dropped from, the equations .further simplifying them U 
Sim.ilarly, and to .first order only, rolling motions will 
not produce any vertical force or, in other words: 

( 2 . 11 ) 



^ 7 . 

6 



6Z 


6Z 

- 


5£ 




6K 


6p 

0 


0 


6v' 


6v; 

0 


6 z ^ 

0 



= 0 
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Also the f^quare termc 6^Z and 6^K will be nonexistent, 

Instead, cons id era,t ion rust be given to frictional and 
vortical effects v;hich are of the form c^^w |w| and c p jpj . 
Although, these terms are not analytic, they will be 
included so that the model represents a closer approxi- 
mation to the vehicle physical behaviour. 



The form of the equations will then be; 

r. m 

X w-r^ .w|w|-t-z, z i r. '■ 

, r V' v/iv.'l ' • 2 n •— v r v r > 

w+- 

!- input 



,n , ,m 



w_ 1 Z w+Z .w|w|+Z_ z +...+ 1 Z„ x'- (x.) (x )■ + 

pi«i ’'o ° rmpp ='10 1 0 

( 2 . 12 ) 



2 - W C 

O 



ios cj) 



-K . 4-X 

P ■ XX 

+in"ut 



K n+K pin] r 
c n I rl M l 



t • t t 



n 

1. _ K n n ( - . ) ( ^ . ) 









X 

2 i 



4> - 



vfnere 




and K 

X 

i 




T 



In the next chapter each term in the above equations 
which is of first ord'^r, quadratic or at most cubic v;ill 
be discussed in relation to its nature, order of magnitude 
nnd m.ethods of predicting its value. 



In summary. 



from Newton’s equations and fer a case 



of restricted motions, assuming the existence of a Ta^/lor 
series for the force and moment and linear dependence of 



I 
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force and rrior.ent on acceleration it was possible to obtain 
S’ set of equations in a state space ropresente tion to 
which some system identification techniques may be 
applicable . 



k 

i'/ 
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AFTER 3 



f'O ? I ON PA RA NETS RS 

For convenience of presentation, the discussion on. 
the notion parsr.eters will be subdivided into thr®e main 
groups, namely: 

i) notion parameters resultant of differentiation 
of force or moment on the body with respect to position 
(7^,<|)) only (hydrostatic coefficients) 

ii) linear terms resul tant from dif f erentiation of 
force or moment on the body 'vith ?respect to acoeloration 
and velocity (added mass or added moment of inertia and 
damnira coefficients) 

iii) non-linear hydrodynamic terms 
3 . 1 - Hydrostatic Coefficient s 

These coefficients are dependent only on the geometry 
and mass distribution of the body, for examp ?.e, shape 
and relative position of th^ metacor. ter, center of 
buoyancy and center of gravity, for the restoring moments. 

The coefficients are calculated based on the li'^"'S 
plan of the hull form, 

V/ith the ship in the upright position, the huoyant 
force B (z ,0) is calculated for each waterline (each 










..i : 



I 
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seleotod value of z^). Then at each waterline an angle 
of roll is simulated in the body plan, being the point 
of rotation the assumed position of the center of gravity, 
and the resultant buoyant force 3 (z^,(j>) and the righting 
arm GZ{z^,<^) are calculated (for details see, for example, 
"Principles of Naval Architecture" edited by J. P. Comstoc'k, 
Chapter 2). The righting moment ( ) is then 

given by the product of B(z^ 4) • GZ(Zo,<f). The vertical force, 

positive dov.'nwards, is the difference between the 'weight 
of the body and the buoyant force, or in other 'words, 
considering the component of the force along the body 

n y T p ^ YJ n I"' 0. 0 j 



K(Zo.f) 



(V/ - 3( 2 q cos c> 
~B(z J>) ' GZ(y^o>h 



( 3 . 1 ) 



where W is the 'weight of the '/ehicle. 

The values of Z(z^,<j>) and K(z^,(f) are us^jally 
presented tabulated or in graphical form (see Appendix 1), 



In order to get bhe hydrost?,tic coefficients, derivatives 
of (3.1), it is necessary to obtain a mathematical 
representation of those surfaces, for example, a pcl;^cr'.omicl 
fitting by the method of least sc'uares. One such scheme 
'was developed by Kerv/in r^l end implemented by Fai'issis £?1 
using I.egendre PolynonieJ.s to represent ship forms wh^ch 
is equally applicabie to (3-1 ). 
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If o'^. 9 is int orrs ted 5e using in the state eei:ationG 
just the linear terns such a procedure is not necessary. 
It is easily seen that 
6 Z(zo,^) 



6 (j> 



4=0 5 7 . 



= 0 



)-fo 



(3.2) 



due to port-starhoard syrnn.etry of the hull forn and that 
6Z(zQ,if)) 

Zo-4>=0 

(3.3) 



6 z 



-^gA 



h'L 



5 K(z^,(j>) 






(j>- 0 



v.’here A-.,j represents the area of the '.vaterniane a.t the 
s i 'n*i waterline ( 0) and is 'y*e inetace^tric i^eiept 

(distance from the netacenter to the center of gi-avity). 
Due also to port and starboard synnetry . in the non- 



linear formulation the term 1_ 6 K(z , ) 

2 ! ' 



6 Z 



IS 



4^0 



equal to zero. Therefore in this case the cubic te'cn 
'■•'ill be maintained in the '*qua.tions. 

3.2 - Added f'asst Add'i'd bonent of Tn.'^rtia dnd Da'^.nina 

Coef f i c ient s 

From the sight tei'ms that form this subgroup four 
’.'.’ill be zero again duo to symmetry of the ve}iie''e. Those 
v/ill be the cross c curling terms 



-Z,.=z _-.-K.--K -0 

p p V.' w 

and the terms to be r^tam._nod ; 



(3.-’t) 
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-Z.' 



_7 



W 



-Kj 



-K 



■D 



( 3 . 5 ) 



So it is seen that v«'hen only linear terms are 
considered the equations of heave and roll are completely 
decoupled and may be studied independently. If the six 
degrees of freedom were considered, the existence of 
port and starboard symmetry would cause, v;hen retaining 
only the linear terns in the equations of motion, the 
decoupling of the eo^uations into two groups 

coupled surge (translation along axis), ■'•eave 
(translation along the x.^) and pitch (rotation about the 
y.p axis) - the so-called motions in the vertical plane, 
and coupled sv/ay (translation along the Xp axis), roll 
(rotation about the x-]_ axis) and yaw (rotation a.bout 
■^he axis) in the hcrisontal plane. 

Added m.ass and da.mping coefficients are in renera.1 
not known except for some, special types of idealized 
forms and when the vehicles are excited into some special 
types of motions. One such type that have been studied 
quite extensively, especially in the domain of surface 
vehicles or vehicles near t'^e surface, i.s that of a 





# 
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vohicl? oycited by vjdv^c into an onoillatory notion, Tho 
path nsnally followoii r3j'°i°1 i?" to anaurne the fl’iid to 
be inviscid, inconpressible and irrota tional , therefore 
assuring the existence of a velocity potential that 
satisfies both the Laplace eouation and the Bernoulli 



equation. Considering the notion harmonic, the solution 
is sought as the product Re L'f ( ^ ® u; is 

the freo_^uency of oscillation. Considering the linearized 



problem, ^(x^,X 2 »x^) is assumed to be the sun of eight 
other potentials being the total notion the superposition 
of each due to each of the potentials: t'vo of them due 

o i 0^ Q 1 0 c ^ v/!?^ VP p 1 u P d i p? p 

on a 'Tixed vehicle and six due to forced oscil la.t ions 
of the vehicle on calm '.vater. The fluid is assumed in 
all cases to extend to Jno in the xq and xg directions and 
bst'.veen x^-0 and or bet'vcen x^~0 and x^=h for the 

finite depth case. The linearized bo-^'^dar^'’ conditions 
will be: 



a) On the surface Xo~0, the velocity of the fluid 
wheun 5 is the surface '=1 eva tion , all the 

6x^ 6t 

potentials and is equal to 'P- using the lir'^arix'^d 

p* ' 



Bernoulli equation. 

b) On the body for the incident and diffracted 




# 
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rotent ialr = 0 = v/h pre n is the unit norpnl, 
6n 5n 

For the six other potentials th';^ ’^or^al coinpon'^nt 
of the velocity will enpjal the velocity of th'^ body in 
that direction or _to+ . or being 



s . =/5 - iojs 



• i‘*?t 



.la 





_ • 


5n 


= f -S . 

.1 3 



(s • is the arnnlitude of the oscillation 

'no ■*- 



•la 



in the j direction) 



6n 



= -ioi>f .s . 

T ,ia 



(3.6) 



where f. (j = 1,....6) are the generalized direction cosines 
J 

of the body surface 

fj (3=lf2.3) cos 



f 



4 



4 

c ) At 



= +« 
= +h 



e _v e 

■'2^3 -3 2 

bb-ns 

b-2-Vi 



V<^^- 0 
$ _ 



(3.7) 



5 ; 

6n 
2 1 



0 



d) At r = (x^ +x^ ) 



the radiation condition ri’st 
be satisfied for the six potentials and for diffraction 
potential or: 



lirn I’" 
r-*-» 






6r 






'■r fs 



= 0 



j = 1. . .6,d (3.R) 
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Alf^n from Eernoulll eQ.uo.tion the hydrodynamic pressure 
V.’ ill he 



t) = - 



r 



6 

6t 



^rr Pp ' ^ ^ = 






(jO 









( 3 . 9 ) 



The hydrodioiamic force in the k direction v/ill he 



= -ff pVs=Re 



:rv/Tri,fdv^ 



s- 

0 o 

The total hydrody'"'amic force on the body in. a given 
dir-'^ctien k was de'^ined to he in Chanter 2 for motions . 



( 3 . 10 ) 



in mode j: F, . - 6 F^- s. + SFv.; 

" kj k^ .1 kj 



J 



s . 

J 






6 s . 
1 









3 



rrz P r:s r.- ^ 



• - <A>t - 



9 



s- = R"^r-o's. e 
J ~ 3 ^ 



--i‘^t 



and changirntT the nomenclature to 
6 F 



= - 



6 s . 
1 






— i_ii. = _ u 



6 s . 

O 



kl 



P _ 1 ^ 0^13 

F, . = Rer(^ A, . + i w B, .)a. e 
k J - k ,T kj js. 

an^ the total force in the k direction vas 

6 

F, = F,. . + F, 



( 3 . 11 ) 



k 



3=1 



kw 



F^,„ the excitation in the k direction (the hydro- 

static force "’as 'Excluded from the .foremoi P" OXTT''^ 0 G s ^ 
as it was from the Bo’cnoulli equation ( 3 - 9 ) ). 
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Conparing the foregoing exproseion with (3»i0), the 
f ollov.’ing is obtained; 

( 3 . 12 ) 



_ TPo 

^v/k - •'- 



(<0A,. + ’ -J-'fXi' f- 






( 3 . 13 ) 



and therefore v:e obtain 

u> S 



■V: = " I" H 









^ - P /J 

s • ' S 



( 3 . 1^ ) 



ja 0 

that would give the value of the hydrodyr.'=rnic coefficients 
:f th'=‘ velocity rotentia.1 was kncwn. Since that is not 
the case, it is n'ct possible to use the a'''o\"‘^ ecuatirns. 
Instead one looks for possible relationships batv.’oen the 
hyd rodynar.ic coefficients to get a better insight on the 
solution of the problem. 

For example usi.ng the second form of Green's theorem 
to a velocity potential end to another velocity 
potential with k 7^ j v.’c obtn,in since both pot-^ntials 
satisfy Laplace’s equation: ( 3 ! 



A,,. - A. 



k 

B, . = B, 

]<: 1 l 






( 3 .^ 5 ) 



Al^o Tron [^^"] n.nd [*10 j thn rolldv/ln,^ ox'nrop:o .Ion is 
obt.oin^d r^^latinfz ^dded rrsss o.nd donoir^?: ooeffic loots : 
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oO 






TT o_ 

OO 



t - 



(3.16) 



- 3,.(H =|f,[A,.(t)-A,3(-)]i!.^ 



r'layte the nost importa'^’t relations are the so-called 
Haskind relations [11] based also on the second form of 
Green's theorem when applied to equation (3.13) after the 



expression for f^ in (3.6) or: 



e - 

■* Ir 



'.lS 

in 






ka 



is substituted in (3.12): 

- ^ cot 



■p 



MG 



, v' o n 



v/k 

ka 

’ tb on er' 



-f- 



r 



f or 



fa 



^ d: 






s ^ 



( -r'f^ ) 6^< ds 

6n 



(3.6) 



(3.12) 



' 6n 



H- i_?) 

6n 



(3.13) 



Since both potentials satisfy Laplace's equation the 
loft hand side is sc^'o and therefore! 



Sf 

’ o 



vp 6lk 

b' ^ 



ds - ff -P 
s ' ^ 



i]5. ds 

f-n 



k 



(3 , ] A) 



ind from the boundary condition on the body 

6*]^, = - efv,’ 



(3.15) 



6n 6n 

The follo'viny will result for the excitation ■f'orce: 




1 
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F =F‘- 
^wk •' 



45_ fS [f-; ^ A] 

" 6n 6n 



L"ka - o 

olirriina tin>^ tho diffraction -Dotontial. 



(3.16) 



Follov.’ing Clll lor 2-d inons ional bodies (under beam 
seas) the following relations are obtained, based on 
(3.16): 



Pi 1 - Ul> 
kk - — 7-2^ 

and therefore 



ka 



(3.17) 



Fv= = C_£3kv] 



r:a 



= 2,3 J 



u> 



where is the amplitude of the incoming ’.-.’aves and 
is the amplitude of the exciting force. Being the 

^ 3'' ^ i R is 

at in'^j.nHy it's nc^'sible to calculate the ma-rnitule 



of +he exciting force. Hov/ever, it's phase is not 
determined. 

From the foregoing it is seen that the analytic 
solution for the genera.l case of a three 'dimensional body 
does not exist. Only the cases of the spke?’es and 
ellipsoids have been solved C12J. H ’-vould be possible 
to solve a. problem nume.rica.lly by such a procedure as t>e 
finite element technique but nothing is av'^ilable yet 

In order to get information .for practical application 
3 . t is necessary then to v.'ork into two d i mens ion.al prob.lems 
The strip theory for ship motion.''', introdu.ced by Korvin-- 
Kroukovsky Cl3» 1^'m -5~\ makes use of two dimensional 
.aol3;tions for each transverse section of the body 



d et-^r’ilning the characteristice on the fl.cv; end conrre- 
quently the local hydrod^Tiamic coefficientc and then the 
res') Its are integrated over the length of the "body to 
obtain the hydrod;y^arT!ic coefficients for the whole 
vehicle, applying then a correction factor to take into 
account the existence of a. three dimensional flow. 

The solution in two dimensions has been obtained 
for the three possible modes of motion: sv/sy, heave and 
roll. One v/ay to obtain the solution is, '.jsing conformal 
mapping to transform the cross section of arbitrary shape 
cf a t'.7o dimensional cylinder into the unit circle, 
to as'"”^9 the ■'/elccity potential to be ^ sum ef ’■'■u'' tipole 
po t *^nt T a.ls whose coef f ^ V ents arc adj'isted to satisfy 
the boundary conditions on the cylinder For 

example for heave the multipole expansion '/.'ill be formed 
by the sun of a source potential at the origin plus the 
m.ultipole potentials while for sway and roll instead of 
a source there will be a dipole since the heave potential 
is s vTnr^ ctr* icp- 1 *th 0 o'^r-hGi' 'tv.^o 3-svTi'^G'tr‘icG.J. * 

The conformal transf or’i’S t ion will be of the -^‘orm 




N 

E 



(2n+l) 





(3. IB) 
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which maps an arbitrary section in the z plane into the 
unit circle in the ^p]an<^. From the section offsets it 
is necessary to find the coefficients A method 

of determining; those coefficients ’.‘.’as developed by 
Bermejo [l?l using Theodorsen's method. 

Another way of solving the two dimensional problem 
is to represent the velocity potential by a distribution 
of sources on the contovjr of the section [l?!. T’^.e 

intensity of the sources is determined applying the 
kinematic boundary condition to the underv'ater portion 
of the section. Once the potential is obtained the 
nr'^ssure is calculated us^ng ecua,tion “ linearized 

Bernoulli equation - and the force is obtained by 
integration from where the added mass and damping 
coefficient are calculated (Eq. 3.9 to 3.12). 

Based on this method a computer program was 
developed []l9l which perform.s the above ca.lculat ions 
and 'vas used in this investigation for the computation 
of the initial conditions of those motion parameters. 
Listings of those programs are presented in Apnerbix ?. 

Still another way of calculating !-!y Irodynamic 
coefficients is to use approximate formulas, based 
some of them on the same principles as 



exposed above, 






I 

?;■ 

r 




» 
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and some in model experiments and full scale infermation. 
Some of those formulas are presented hy Blagoveshchensky 
[ 20 ] and were used here to calculate » ^33 

®44 results compared v/ith those obtained w'ith 

the computer program mentioned , believed to be more exact. 
Those results were also compared v;ith the ones obtained by 
Hines and Gies [2l1 in experiments v/ith another model 
of the sane hull form , Those results are presented in 
Appendix 3. 

3.3 ~ Non-linear Hydrod;'n^amic Coeffici‘^nts 

In this category are the terms of second order in 
velocity and cross t'^rms in. velocity a.nd 000 it ion. 



In.f orm.at ion. on those types of parameters is 
p?’a.ct ical ly nonexistent. Only two references to 
frictional damping in roll were found [20, 22^. Si.nce 
that is the case where frictiona.l and separation 
resista.nce ha.s a bigger effect when compared with t]ie 
other modes of motion, the expressions suggested by 
Blagoveshchensky was used since it ’.vas directly applicable 
to ocean vehicles. 

The expression given by reference [22^ relates the 
friction coefficient of oscillating cyclind.ers to an 

f 

equivalent Reynolds number: 



= 1.328 



'<.22 9m 



Tt) 



-.5 



(3.19) 




# 
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v/here y- ~ radius of the cylinder 

9 - mean anrlitude of the oscillation 

m 

T - period of oscillation 
- k.inematic viscosity 

This formula is the sane as olasius solution for 
in laminar flow if the Reynolds number is taVen as 
Rey = (3-20) 

To obtain the drag coefficient foi' heaving 
oscillations a sim.ilar m.ethod may be used, although with 
rougher approximations since the vehicle is not a two 
dimensional body and the sections do not resemble circles. 

o ishi ri'P'D^ro'”' c^i v/?s R0*nij'Li3s 

were sought fc^ each station and integrated over the 
length of the body. 

For each station it is assumed that the resistance 
is due, when moving into the water, entirely to skin 
f'ciotion and when moving av/ay from the water, entirely 
to separation. Therefore for the first ease the 
diameter of the eogi iva,lent cylinder will be d = i/r 
whnre 1 is the length of the underv.'ater nart of the 
contour of the s'^otion. For the second case t’-^e diameter 
of the equivalent cylinder will be equal to the beam ^ 
of the section. The velocity used, to cclcula.ta the 
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Reynolds number '.vill be the "root mean square" value 
based on the tests performed. Values of skin-friction 
coefficient and form-dra .3 coefficients are obtained 
from reference C23D* Results of this procedure together 
with the results of the coefficients as calculated by 
the method indicated in reference [20I are presented in 
Appendix 3- 

Vi'ith respect to the other coefficients, the 
assumption nay be made that they are indeed very small 
and consequently the assigned value to bo used in the 
identification procedure is zero. 

S u IT ^ i *n ^ '1 3 "P 

\v0]r0 G 3 03rGd > i'h.GiT*^ nn. t'' • ifG j ^ Tr.'"oo'*'''*^*0.3^GP '^nd 
the available methods of prediction outlined. The results 
obtained for the chosen hull form, are presen+cd in 
.Appendix 1, 2 and 3* 

In the next chapter the exciting force and moment 
•will be considered thus completing desc:>^iption of 

th e ma th ema t i c a 1 r. od e 1 . 
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EXCITATION FRO^' T':R S:' VTRONITS?’T 

Froir; all pousiEle tyros of oxoitatior, o" 
and ncnents duo to tho action of two dinens’on 
progressive wa.ves are considered lecau.se (l) t 
v/ell defined in mathematical terms and (R) eas 
generate in a to’wing tanF experiment. 

The exrression for- the ve'^ocity potential 
obtained following a rath isimilar to the one f 
]-v ^ p ^ f'ojT c P Cl 1-3- i”/ i on oi* 

linear hydrodynamic coefficients, naTT^oly tha.t 
"’^locity r o t e*'^. t ia*^ 3^tiefi' ?e hoth Lan^’^ce an^ 
eo^’.’O tions . Takvna the case of an infinite flu 
directions the fol].cv’in“ exprcss'’ons hold ^ 15 » 

,, - 0 



d(f-^x )-C 






V J 



aM-m 



0— (4-,- ... -1 -^5) 

ft 



= 1 ifc: 

r 6t 






X -0 
3 



ly forcec 
p 1 r lac ^ 

V' p > r p » Q 

. j 

y to 
is 

ollcv/ed in 
7 /een the 
the 

Qi ; *1 1 ^ 

id in all 

oa1. 
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y — oo 






= finite ns r =-/ t- x' ^ 

JL 2 

V 'h 02? 0 is the surface displacement. 

The solution is sought by sepa.ration of variables 
and a.ssurriing a harmonic dependence on time, or in other 



words, and for two dimensional behaviour; 

-iti>t-| 

) I 



(> 



.^2 .>^3 .t) = -e[^2^^2' 



(^.2) 



Plugging this solution ba.ck into Laplace's equation 
v/e obtain 

^"(>' 2 ) = “^"(^ 3 ) = ^^-3) 

~ Y 






3 ^ 



To satisfy the last boundary condition ^ust be 
positive being the solution for 

ikx. 



fd^> = e' 



■2 -li<X2 

+ C e 

2 



kx 

D e ^ 
1 



+ D e 
2 



^(x ) 

|3 3 

?nd to sa tisfy y <p = 0 x — ' D 

Tw 3 1 






= 0 



f (x ) 

|3 3 



and thornfore 



D G 
2 



-V X ' 



-kx. 



S (x ,x ,t)-Re [ e ' """(C. 

I\\r 2 3 



i (kxs- «'^t) -i(kx2-:-«Jt) 



-C G 
2 



) 1 ('•■.'5) 
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Th"" '■olut‘''^n v.'ith C _ 0 corrocponds to '.vavGc: 

2 

■oro?rresG :i np in the oositivo >: direction and ”,-ith - 0 

3 1 

to waves prorroosing in the negative X2 direction, 

Assuning the later cane and the existence of an 
erhitro.ry phcee angle ^ to he determined la.ter from 



the initial conditions; 

/Tv / o /, -[kx + i(kx f uJt +6 )“| 

Q) (x ,x ,t)-Re C e 3 2 -* 

J W 3 3 



(4.6) 



lomhining the tv/o houndary conditions on the surface 



( X 0 ) 

3 



I 

c 



6t 



5 V/ 
6x 



X = 0 



CO 



= k 



(4.7) 



'vhich determines the value of k, the 'vave numher. 

hsing 1 o$x 

J " g 6t 



-i(kX2t‘^4 t 6 ) 

1 Re r -iOJ C e ^ 1 



wd sin (kx + >ot-f6} 
2 



h . 1 



bvA f or th 0 f or^n 

*rr sin(^x 

'v'!"'. ^rc ^ ^0 1'^*'"' o or.'pJ 3 *^■’7 do of tlio our^o.C'^ oplp. c 

,n ^ O ^ 



n 



CjJ 
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anc 



$ (x ,x ,t)^Ro 

' 2 3 



— rv-.’ -L. * co 5 ) 

e -‘ 3 * - -^'2 ^ ' 



tO 



(/i. .^o) 



from the Hashind rolatiorr, So. (3.1^) 

f, - -f? ] -I- 



P — Oc 



C! 



s 



6n 



6n 



\V o 0 



'Pb (O f h *t ■^'GX'P T p "t "1 “(*V> r '^or'0 0 

exerted "by the pressure in the wae/e v;ithout teh^rg account 
of the dicturtance by the presence of the ship o.'^d the 
ontv excitation accounted ^or b^.'' '^he Frcudo-Kr5.] of f 
yoo th es is * The seconxd tern is the ccntr *hotion of the 

T p C '■^'^ O t ^ ^n C t ^ \\tp ^r <r* py^ OU '^b iOl r> b r ^ 

Tc oormute ^-ho f^^rst terT^ in the n^c n 
(the Freud o-Kriloff force) first replace 



t . 



1 tJf 



6n 






(3.6) 



V.'.l fK 






r r 






n 

^rr e' 
S 

o 



[hx^-^ i ( Vx^-{- t-^6 ) 1" 

f . d s 

1 

o 

3 ^ T n. ( hx -f 5 ) f . d s 



e' 

lO 



(^-•'.11) 



The secord term (cl iffractior. “^orce) i r ccmpi’t 
notin^^ that 



b% = 
6n 



, I V y ^ ^ ^ h X 6 ) 1 / . X 

3 ' 2 (f-f-if) 



( T? ) 






■f— 

3 a 



- ^ OJ “ 



f. 



• ft >y 4-"^ ( h"'" ”^*6 ) 



(f +if )dr 

2 3 



b' 



(li 19) 
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an-i 



v:i fk v;j'd '• * * 

Th^rGiorc ‘to oo-lculs-to "the tot^l oxc i'tsi't ion rcrco‘ 

the results fro"n the eut-prohlem of ‘‘'"orceo cso illatior^s 

of the body in calm woter must be kno'>vn , 

4^(x-|,x ,x }. V/hat was said abon t t'^e v/ay of calcu ] ing 
added me.ss and darrioing coefficients .is npo"* icab3o hena# 
Besides it is Possible to calcu.3.ate this second tenn in 
tt e exciting tQ 2 r‘ce at tl'ie Sc.nc time as the edded ^'^pss 
and damping coefficients ane ca.Icnl2ted i since the 
potential is also obtained in th^ node of interest, Cr*, 

^ ot * "\'0 obt'^T'^*> fo"^ o4-p-i-TQv^ 

■ CJ - [_ky 3+ ■’. ( ;<x 3+ 5 ) "1 






9 -^ 1 >71 






( f -r i f d 1 



{h . 15 ) 



v/here C is the contonr of the station and d] ^s'the 
o 

length element along th'^t contemn. 

These two dimensional results for each s tat*' on ore 
-♦-'hor g“'^a t d o v c> 1 opi rr*[' ”d o ^ f' ^ 1 e , 

A 3. t h on a'o t h i s e t o d s e m n t o h t h c c n t s ■* ' i t e d 
for the commtat^on of thp pxai*^^na fc^.'ce^' 1 t^’^'^t w'^s ^*'c 
u^ed in this invest ira*^* ion due to t’^e '"'hat '^he 

computation of eaination (^’■• 15 ) '-vrs not i^^.c"^ ".d "^d ~’n *th 
0 omen ter scheme used ir the ca cu. la. 1 5 on of sGctiona.3 
a.dded mass and dai'^ping coef f ic ie^'its . Ii^s toad , the 



» 



'M 

1 

'i 



c 
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nroceivire precc-nted "by BlJ’-y^vorhc^enoV'y ^20^ v.'as used 
since, elthouph less correct, rendered the cor''puteticns 
ruch easier to perform. 

First it is assumed that «1, , ?/;;\ «1 

•X 



■/.’here 5’i is "^he v/ave amplitude, ^ is the v.'a.''."e length 
A“ Sm/h , B and T are the beam and the dra'f't of the 



vehicle . 

Next, ^t is assumed that the principal part of the 
'^xcitetio" is due to change in bn.oyancy due to the 
presence of th‘^ Vi'av'^, Then some multiplicative correction 
coefficients a.re ca.lcu].a + ed to account for ■*'h'^ curvature 
of v/ave (finite hear ccrrection) and for the den‘~ndence 
0 "^ t ^''0 h’^d ■rod'''T.c.mic 'oressure on x (tir>ite draft 



correction). Fach correction is calculated, separately 
and assumed independent from, each other (for exam.plf^, 
v/hen calculating the beam correction, T,/^«l end vice 
versa ) . 



Under these assumptions the ’’ inear efunt.^ons for 
hoa.ve and roll v/ill tabs the .term: 
a) heave 

V/ — 7 \»r \- 

V! * 

-Z *+111 -Z *trn , , , ^ , 

W V/ (^i.l6) 






V/ 







A 
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If]'' o V* 



fc ic 



njr.f'?<C 0 c is pl 3 .c ccr'r'ost^i Tex' 



coan ari '^ra.ft ef.'^ects 



1 = 1 - c 



b) roll 



P = 









-K’+I 

p XX 



zb 

sin (ky + t + 6 ) 



"Tir-'+i 

V XX 






4 



(^^. 1 ?) 



-K- + I 

p XX 



v’hsT '0 is 'fas v/avs slst)0 coxi'oc't.Gd fer tsa.ri and 

d xa f t ^ P ^ p c t s 

-f - i.. ■ " 



Vb 



Vt 






cos(-^y ^ i 'i -6 ) 
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The detpilc of t?.e calcula+i ors of thor:o corroct.ion 
cooff icionts nay be found in Reference C?o] para,<^ropbs 
27 to 31. 



the contribution due to the inclination of the chip v;hen 



will not be in the same vertical thus causing a heeling 
moment V/.EG sinc|> being V/.EC-<|> the first term, in the 
Taylor’s expansion. W is the weight of the vehicle and 
,EC is the distance from the center of gravi.t'''' to the 
center of buoyancy. 

,A.lt’''ough the '^xoress ions a"^ove for the exciti.'^a force 
atni moment ere not as correct as Equations ^.1]., ^J-.13 
and 4.14, the simplicity of their calculation decided 
their introduction in the equations of motion. 

'A’hen considering the non-linearities, if the 
exciting force can still be app'^ox irnn.tcd as expressed in 
Squ.ation 4.l6 and 4.17, the follov/ing terms should bo 
a^ded to each pair of equo.tions. 
a) Heave 



In Equation 4.17 the term V/ . Eg 



rear es ents 




ince the weight of the vehicle and its buoyancy 




V w 




-Z • + m 
w 



where the first term rerr'^sents the tv’iop? or.ol and ■f'orm 



dra.m and the other term is the second tern in the 
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polynomial fit of tho Z 




b) roll 



1/2 . P|F| 



1/3! fk.3 + W.EG) 



^ and 





(4,20) 




V/here the first term is again the contribotion of 
frictional and form dreg, the second term is the 
contribution from the polynomial fit of the curves of 
stebil.ity and the last is the following term in the 



for the excitation is even more Questionable d^ue to 
the small orders of magnitude involved in the identifi- 
cation. One v;ay of obtaining a better approximation 
is to consider the effects of the presence of the wave 
in the da.mping terms a.s Is done for the case of roll 
by Tamya []2Z], 

Those terms would then look l^ke 
a ) heave 






(4.21) 



\ 









% 
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b) roll 



V 



V 2 k ^ (p- fO|p- (4 



-k' + I 

P XX 



-k • + ! 

O XX 



wh ere 



6t 



fv = 



"Vi 

6t 



V/hon considering the non-linear coupled h''a.ve and 
roll rode, all the coupling ter"s as nentioned in 
Chanter 3 bad to be added frcn v/hich only the curs 
hydrostatic terr- could he ry'^dicted by using, for example 
the method p'ccscnted in ?.e.'^^r'^rccs^''e , si, sc’^c 

cbyections v;ould .apply with r''sp>'ct to f' c>:cita.ticn. 

In summary, the excitation due tc the ^.ct’o:" of 
plane progt'essive ’.vaves can be calculstcd cuit,-^ accurately 
using the Ha.skin.d relations. Icing the action, from the 
diffraction potential calcula.ted for each ciatio'" 

"tv/ 0 d. in m fa 1 0 1 i- 0 c i n \ ^ u 3 s -^r d d ' 0 n 3 3 u 1 ^ 3 n "t ^ 3 g 3 nd 

nonGni; o^oi;.3.ir;"3 vv int0Tn3.*tion of fdoso di**o ^ ^ 

]r 0 3 n i ’o s 3 .1 0 n n t :i 0 Jon n *' ''i o ]i 3 v 0 1 c *1 0 . Dug ~t o 

V ^ r-^ p ■}“ V r' rr p c- o « y'r' ^ p V 

0 0 0 0 0 C U VP- f 3 Vu f off 0 n3 

V 

i t\y Another advantag'^ ■^ot 
n 30 . 11 . i on 0d f of 000 is fu.s/i f ^ s "^otoic '^■'"olrl’iC'd "‘xm'GSs 



compu 


tati ona.1 


d if f iculti^"^ 


(p* ' r p ^ 


■K' n ^ 1 ■! 1 “^0 
0- J U — ^ -J 


believed tc 


the a 


d VP ^ p ^ 


of i":s si-'-.nl 
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the exciting forces and inorrents in terms of instantaneous 
values of surface displacement and wave slope which are 
easy to obtain during the tank experiments setup to 
measure displacements and velocities that will be used 



in the identification procedure 
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CHAPTER ^ 



EX PER I MEHTA L PROCEDURE 

A model of the "mariner" type hull form,, whose 
characteristics are giver in Table 5*1> ''’'S-s tested in 
regular waves measuring the displacements and velocities 
in heave, roll and coupled heav*^ and roll, and also 
the surface eleva.tion, with the purpose of identifying 
the hydrod;^ma.mic coefficients in those m^odes of motion. 



TABLE ?,1 



MODEL C H ARAC TERI S T IC S 



sc3.1e (model/ full size) 
length (L) 
team (B) 
draft ( T ) 



1/96 
5.50 ft. 
.79 ft. 
.28 ft. 



weight (W^) ' 4 - 6.05 lb- 

mass (m) I lb,sec^/ft. 

The model was placed in the miid-length of the tank 
with its centerline perpend icui.ar to the centori.ine 
of the tank such that, using the coordinate ■^\''stem as 
defined in Claapter 2, the direction of p^'opagation of 
the waves generated at the wavemaker in the hegirr.ing 
of the tank, coincide with the negative x - y direction 



The model v/as ola.ced at the mid-lenath of the tank 
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as a sol’jtion of corprornise : not too close to the 
v'avenakor in order to avoid the action of the diffracted 
waves fron the model upon reflexion on the wavemaker , 
not too far from the v/avemaker in order to avoid 
interference from the walls of the tank. Reflexion on 
the opposite end of the wavemaker was of much less 
concern due to the existence of a "beach" composed of 
stainless steel shavings with a reflexion coefficient 
of about . 05 > for normal v/ave lengths. 

Tests were run at three different frequencies. 

First it was decided that those frequencies v/ould be 
half, tv/ice and the natural undamped frequency for both, 
h G V -? G n rt T o 1 1 j v/h g o “t h g g g. t/U vs. 1 xind. g g d cej \x g r c 

is defined as for a linear system 



CO k/m 

n 


where k = ~Z„ - - 

zo 






TTi ^ -z;+m - 





For roll the measu\’'ed natural period (from free 
oscillations of the model in ca.lm water) v/as T “ 1.52 
seconds or f„ “ cps , Therefore tests were run at 

= .33 ops sT^/Rtt = ’^7.36 ft. 

f “ .'^6 cos = 1 1 . P 4 ft. 

2 

f 

3 



1.32 cos 



2,9^ ft. 
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For heave the neaeured natural period was - ,7? 
seconds or = 1.28 ops. For the tests to he run at 
twice the natura.l frequency f - 2.hP cps, it would he 
too close to the maxinur. at 2.50 cps pernitted hy the 
wave rr)a,her. Therefore, it v/as decided to run the tests 
in heave at one fourth, half and at the natural frequency. 



or 

f^ - .32 cps A = hQ.QZ ft. 

f = . 6h cos = 12 . hp ft. 

2 

f = 1 .28 cps - 3.12 ft. 

3 

For the tests in coupled heave and roll, it was 
decide^ to run at the nat’;ral frequency iii heave ai^d 
at half and at the na.tural frequency in roll, therefore. 



for the coupled )node 

f^ - .33 cps 

f = .66 cos 

2 

f =1.28 cps 

3 



\ = ^7.36 ft. 
A = 11.84 f+. 

A = 3,12 ft. 



The quantities r.easiired in each tost •vere 
a) heave 

surface disolacenent 5(t) 

heave displacement ^o(l) 

heave acceleration z (t) 

o 



-51 



The ■^'?3son for' o oo c*^ 

velocity .’/as due to the ivposs t hi 1 i.-^.y of ohtainio/T e 
translational velocity transducer, and to not heinv 
able to build a reliable electronic interratcr '’/ith the 
material available 



b) roll 



surface displacement 




roll displacement 




roll velocity 


P( t ) 


Amain '’/ave slope should 


be measured in this mode 


instead of surface elevation 


but the probe used could 






c) coupled h^avo and roll 




s’;rface d isplace^'ont 


^(t) 


heave displacement 




h^-ave acceleration 


zyt) 


roll .d isplacenent 


4? ( t) 


roll velocity 


r(t) 


T]^e sa’”.e otservations as 


HsTd 0 ^ n. ) 1') '* Vcvp 


a ml 

1 

Description of probe.- 't;. 


' ■ I V ' i ' ^ 1’ 0 n r* u p 



p i mul tan'^cus ly surface displacement e’^d •.■•av-a slope 
is available in Reference 
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SurfacG displacement was mes.sna'ed with stance 

type ’vave probe formed by tv;o parallel stainless '^t'^el 
wires of l/8" (iia,neter. Sirce it was necessary to 
obtain' the wave d3 spla.cenent at the veh-icle, the ideal 
position for the wave probe 'would be along; the renter] ine 
of the nodel. Since at this position the probe would 
be subject to interference between the n!odel and the 
tank 'walls, it 'was decided to put the probe in integer 
number of wave lengths ah-^ad of the model. Those 
distances were calcula.ted based on the relation 
> - o)^ W or gT^/8^ and this procedure introduced 
n firsh c^ror on bhc T;op ‘^ircc f .r ^*nihp 
d>or*ihs rdc.hion b^iv/cp”^ rc'^'‘'hr'T* 

frequency is [23l, 

(0^ V t^rh(i-'h) (5.1) 

T.r' 'o-pw'i "thn 5vs.tor dnnbh. 

For ey.amp] ^ , considering h ~ ft, 

f = 

,5 aps .098 

,65 cps .025 

7.^5 cps ^ 0 

From the above it is seen ‘'or the low 

frequencies f< ,5 cps the er'ror heco’n^s apprec i-’ ble , 



( 



f 



m 
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He^vo d ,i nplacencnt v.'as moacured with a linear variable 
differential transforner (LVDT) connected to the !:eave rod 
assembly. The heave rod slides on a sleeve provided 
with ball bearings to reduce friction. Even though, the 
friction was found to be quite high (comparing with 
experiments had with a heave rod assembly provided 
with air bearings) and therefore there will be some 
error involved in the identification of the heave 
damping parameters. 

Heave a.cc elera.t ion was '■"e'^sured with a.n accelero- 
j’neter placed in the m.odel as close as possible to the 
ex?s cf ^:^etation ■^n order to m.inimize interference of 
ro.ll in the measurement of heave acccl crat ^ on , 

Roll displacement '•/as measured by a rota.tional 
va.ria.ble differential transformer a.ttached to a roll 
hearing. 

Roll velocity v/as m.eas/’reci b'.'' 7*ieens oT a rate 
wro , bo'^rowed from the Ch.arlas Drap'^n Labcratoryi v/hich 
1'^'^ '^c^d ^ T'’ ^cdol ''’1 cti 

p^.rallel to the coordinate axeo located on the : edel, 

Thir "^'ete i uned in the De'^e Snhr^.err'^nce S\'st^n:i 

P‘^ 0 joe'll 've f- ah.le o.f nioe in^ enp'nlar ratee in roll, 
niteh and ya^v simultar^'^ouGly , up to ra^tes cf ^'0 depro'^r/ 



second , 



All the transducers 'Tior.td oned ahove, with the 
exception of the rate gyro, are excited 'ey an A.C. voltage 
of 5 volts at a frequency of 2400 cps provided by 
cP^y-ier preamplifiers Sa.nborn model P205A which also 
received the electrical response from the transducers, 
generating a D.C. signal which may be fed to a chart or 
taoe recorder or both. The a.mpl if isi’s also provioe 
means of amplifying or attenuating the electrical signal 
according to its leveT since both recorders work 'cetween 
i 2.5 volts. 

Since the excita,tion necessary for the rate gyro 
i' V A.C. (f “ 4oo cps), the Sanborn or'^ampl Ific^s 
cculd not be used. Instead the following circuit was 
used (nee Fig. 5.1), since a D.C. power supply was 
available up to three amps; 




The 26 V 400 hs excitation ’.va.s '^rd to I ho rate 
gyro. The output from the rate gyro (400 hz, phase 
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,-ev^rRin.n-) was th-^n fed to a d'^noduletor (PRD type ROR) 
to ottaio the D.C. signal proportional to the arguj.ar 
velocity. This sif^nal was attenuated in. a Sanborn D.C. 
anplifier before going into the chart and tape reco’^der 
(see Fig . 5 • 2 ) . 







-400Hz- 

FROM 

GYRO 










4 




1 

II5V 
|400 Hz 


DEMOD. 

PDR 80S 
1 6 


-L.5>^F nr 

-j- ^ AMR 







Firure 5. 2 



For convenience of futu.re CcC culat: on , it was decided 
+ ^ rut the center of gravity at the derygn v’aterline 
^'■1' ‘^tabiliti'' calc’,' la. t i one were trade ■‘I'or an a:':"uried 
■ Kc : ^28 - 1 . - center of gravity '’t the v.'o tori ire) . This 

was performed by trial and error by creating a knov.'n 
r.or.ent (with the model in the upright position, a g’ven 
weinh-i: is "oved transversal ]y a certain given dista’''.ee) 

i'^ n* t]t e resulting i. nclinat i.-on 'iVhn.c}'* s'^S' acmp>^rod 
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v/ith tho ].oula.t ed curv^^e 

%^o5.r^h.tc v/oirr' th. on ^ i! 1*/h i 

verified again urfil reaching the value 



of ^ ta h-i 1 i ty , 
p n •'^ t h e T”! ^ 1 e 
er^eified ^n the 



curves . 

Ch- 1 iorat i on 0 “^ all ins trurontat i on ?/ar na'^e in the 
tank except for the rate gyro v'hich v;o3 per'^ormed at the 
Draper La hora t ory , 

The LVDT was first zeroed hv slidina the core into 
the casing until no o^'^tpu.t was s^. ow>^ on th*^ chart rocec''d'^r. 
This noint wa.s narked on the sr^.aft of the c'^'^e toa*ether 
th marks sna c^d 1** on each s i d , mi- p o 

th ^ con^ ot 1 o'*^d her e. •'ti ■»/<*) p p-i-ov-vT^fi o ^ aur* 

T f i'^'^s t^ie na>' r^o f*! ection c^' th^ c^^rt -"tp r» 



noted (1“ ~ 10* mm - 50 attenuati on -tor) , gn-v* 



ohe end 


i’otat :• one], di 


f f er ont ial '^^’ran 


s t or m 


'P 3 './S S 


cip':i].an, Fot' 


h h e c a T c ^ 


‘t TP s f a r e r 


tz''" c tor 


das* d h o ■’'* 


t a. n* ' *^p o s , 


mu pp 1 - _ 



1 ; o 1 ,rr n r c 

‘'cration cf th? ’^rv? ro-'l "’'(i ■-.'■-1'!, "oearing v/as checked later 
v/ith the instru:nontat ton a’coard . 

Pot* t ^ ^ ^ * h*'''p t ^ o ^ ^ o p r> ^ ^ v-* r* *vu s 



ec^ed w ^ tt ^ t i a t^’'^ o ^ * sot t "* 

v/as thc^^ i^mr^ad ^c T-n ^ and " ts 

( *^ * 5 a* ^ 1- P , a f f G c f ^ ”** sen f*' * 



■^u . iic^ 1 ■'/a In.':' 

o ^ 1 o ^ i •' 0 t i ^n ■''■*' *h 






I 
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the pon (25 Tr>rn) . Then the zero snppreesion control on the 
preamplifier 'vas adjusted to bring the output to zero. 

For the ra,te gyro, the whole assembly (rate gyro, 
inverter, transformer and demodulator ) v/a.s placed in a 
rotating table to v/hich a known constant angular velocity 
was given. The output voltage was then Tr^easured with a. 
high precision D.C, voltmeter. It was noted that 'the 
presence of demodulator introduced some hysteresis on 
the output as well as some non-linearity a.bove an angular 
velocity of 20 degrees/second. On the average the orjtput 
signal 'vas of . 13^V/degree/second D.C, while for 
f'l.e rate "■yro alone it was , I'tyv/degree/s'^cond tggs . 
Afterwards in the tank a constant angular velocity was 
simulated by 'reans of a potentiometer, the output voltage 
read with an oscilloscope and the pen deflection on the 
chart recc'f’der noted. 



Th*= 


resu 1 ts 


of the final 


. calibration 


for al.l 


the 


trune 


n*ts Vv'P3T0 


as follows: 












wave 
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fa c tor 


= 50 ) 


roll 


>aop -p 1 p,(T 


0 

0 
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10 


rr.m 


(atten . 


factor 
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{Tyro 
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10 


rr.rn 


(atten . 
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Frorn the above eo"^e other ristahcs can be notch . 

First, full scale (25 corresponds to 2.5 volt". For 

a nodel of 3«375 irches of draft heave d isplacf^m<^nts of 

more than 1 inch ( I /3 T) can hardly be expected. Therefore, 

the output signal on tape v/ill vary betv/een ± 1 volt at 

most. For the acceleration if it is assumed 'that, as for 

]^n'"'ar systems z = u?^z the arnolitude o'^ i’^e acceleration 
' a a 

•'.'ould be for the highest frequency (f = 1,2P cps 

- 2 _ o4 sec.~^) - 5^!- , 6^1- inch/sec.^ or 5<3^ ft, /sec 

since full scale was 1,5 g = '!'^,3 ■^t,.^sec,^ the signal 
o'o ■*r ^-'09 V9iry '^i)9'tv/9Gn i »25 vol'ts , For lov/or 

0 '' rr^r* ' 1 r* q oc* -Mg n H 0 O n *5 a 01 J. ^ ^ ^ Ci-n '* a '{hi n> 

^ o '7 1“ -oov^ .9 r T1 OiS ’■ - '^9 "FI P “b'" ^ ^ h* * 1 fp +; l ^ r: o 

output signal approochod the noiso level. 

Sirn.ila'oly for roll disp]acement end ve'^ocity e/nd for 
the rave probe, the signal level was guito low. 

? e .f ore ru r n in ,r the 1 0 s t s call h r n t i cn s i *'^r p 1 e w e r 
3 op t into the tare recorder in order ‘^'o b'^ able to scel^ 

i5 q correct vabc'^s of d. i cplpccrnn'^s ^^elce i"^' iec 

enen the d 1 p"! t ] zat 1 on boin<-'^ nerf on^'-u. About ^0 , of 

.. . - 'Oi '■ ^ ■ 

o -n 0 w 0 r* p r* 0 o q pH o{d v; i t h t 'c 0 0 0 c o ^ i bra t ion c i <^n els t ^ 
p c 0 r d s p o p H o o V t IPO h /s ec . 

Then tlie nine tef-^br’ ac t, ion ^d bc'^ere wer'”^ "^pn 
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recor^ing ? 5 ’nont 50 cycles for cech run. In 'betv.’cen runs 
enough time would elapse in order for the model and tV ; 
water to calm down (zero initial conditions). On tape 
the interval between runs was about 20 tt . 

The waves were generated by a wave ra>er constituted 
by a hinged paddle act\jated by an electro-hydraulic 
power system, v;hich was controlled ly electrical 
sinosoidal signals generated by a low ■'^reo^ucncy v-ave 
oscillator, through a hydraulic servo valve. 

The data was digitized from the magnetic tape into 
conp'jter cards by means of the EAI ^80 analog- ccmp’uter, 
EAT 69"^’ da.t''" conversion system and '*.h'^ l "’/^0 I^h' divita .1 
comeuter of the Deoartment of !’'‘^cha.n.ical. En.^" reerin'^. 
During digitization some mere def ic ie’^c i es v/ere. noted , 
besides the low level of the signals cn taoe (the EAT 6 P 0 
works within a. voltage lev'll of t 10. 5 ^'olts in normal 
opera.tion; wh‘=>n digitized these s'gnal.s correspc:^d to a 
standard precision real vcelue of ~ .05) . E^rst the 

m.a.s;net^c ta'ce )'‘.ad, on th- nortions rot recorded, a D,C, 
signal 0'^ about vol1,s whoso oi''igi'o couldn’t 

0 i; o V'^ A 01' ■** p 1 T f '■') i T'J/ 1 rc r> a f’ o 0 *70 

•{*)a o h 051^*5^ of' 0 c O'*' ^ 

i n f j 3T 0 d p c i n 0 o p s on t ''*1 o 1 p o s of h 0 n q ^ o rr', p I- 
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Also the runs were too short to he possible to set 
up the vn.lues of the pot entioneters conveniently. These 
potent ioTieters were used to take off any bias that- 
existed in the rrieesurerGnts . 

For pr?rposes of digitization, these circuits v^ere 
huilt in the EAI 680 analog computer (refer to Appendix 
A and Reference TRS"! details in the circuitry); 

1 . - Control of Sarrnl.ine Rate 

A 'A'avotek ll6 VCC oscillator is used to generate 
a square v.-avs of frequency equ.al to the sampling rate. 

This signal is obtained in analog trunks 0 or 3- From 

tho 'trunks 'tho ,^oos into tho input ox ^ 

in tho intex'fnce t'^'ay of th^ vow^s. The output, 

'vhich IS otto.ined in the same interfa.ce tray in. the logic 
ro'.v is fed into a differentiator to obtain a spihe. The 
output of the differentiator is in turn fed into the 
analog to digital converter (ADC) sample control 
terminals on the "^'irot strip, pre'^erably on trunh 
This ^ivnal nay also be fed to a BDC cou.n'^er ir, Colunn 01 
to check the ororetion of the circint, 

2 , -Si mn a 1 p r o c ^ s s i n a . 

Frcn the tape recorder the signals are fed 
into analog t:^u.nks 14 to 20, loc<ated below the Ivavc'tek 



I 
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OGoillator. Frcn tho trun’-'^ 'v\oar' outri;-^. is r,itua.ted in 

+ 3r'OP''^^; ^1. Q r* * *'* i" O 3. ^ C 

v.'ith the output of a potention.etor (to to-ho off any bins 
on the sipna]). The output of the nurmor is ted into a 
D n't on “t i *t or* on'i fnoo' ~bho oo'tor-’t lono t ■^’p iniro ^*o 
(input ro-rkC'd with a 10). Tho pu^rrcoo of tho ootor.tlo- 
netor is to attenuate, if >-,ec<^sso‘>"y , "*V'o siynol v/hich 
is amplified ten tiTnes in a^'.p’^ if ior, Tho signal 

connos non lent to th. c n o a s u od a c o o 1 o na i o n , i n r t no r’ of. 

ro-inn* into a.n anplifioir, is fod into an intognatoi' to 
oonv'^rt id: to vr]oc;5ty, T'^e outnut the nmyg ipieu 

rs *v^ -I- r> o*'^^o "h p h p pV 0 or] 1 n t A ^0 t }‘nUo ppo-j-p/^ 

TV- -r}-i 0 O^ d 't* p -p T p ^ 0 A- p 'I ' ' * S n r] r ~r a f It o o p ><a ^ 

pT>n0 in^o th"^ i^sT)~^nst tnun^'^s fon t^^ O^S 

o n ^ OPT os c i 11 on^'c ope , 

3 • “* Oi splay Circuit 

p'* 7^0 TV' O; *p 7 va p ^ P SqP r'i'y* 0 ^ ' ■'O ^ ^ ^ P >>ep^ ^ 

soritnol tinen (in ooluTsn 0? nf tho lo^i^ no’'v) s "‘.'iro 

i'"'/ connnoted to the ^''^cono" — p^ S "’al oosto''! 

rpor* 

in the finsh onden to oh'^''a^n a *u '^ny 
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ccnruter prograin butlt a-ccc^diny t-o tho hy'rnid oGor*s £'),ide 
pijlDlished "by the Depertr'ent of !''’oc’'^anica'' Ve^^ri nc^ 
using existing routines. '^his ^oxr:"ter progren is 
presented in Appendix . 
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p su^''Tr'^ r\'' th rne *- p "no ^ i ^ tn ^ 
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ncip/ted out, If tiioe 'oorpiihod they pUoT-id hr roneatod in 
ordor to obtain h^'^ttor data, 
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CHAPTER 6 

IDENTIFICATION OF SHIP MOTION MODEL PARAMETERS 

Most of the available techniques of parameter identi 
fication have been developed for linear dynamic systems 
[l]. Some of those techniques, assuming the linearized 
equations of motion to be a valid description of the 
behaviour of an ocean vehicle, have already been applied 
in the determination of ship motion parameters in heave 
and pitch namely, the frequency response 1^271 and the 
step response methods 

In the case of a non-linear model, two approaches 
are open, both yielding approximate results. 

The first, making no assumptions about the character 
istic of the uncertainty, assumes a mathematical model 
for the system and then comparing the output of both 
model and systems to the sam.e input, searches for the 
structure of the model that minimizes a function of the 
error betv/een the two outputs. An example of the fore- 
going is the model reference technique^ For the system: 
Yg('t) = ,t) 

y (t) = f (u(t ,t) ,t) (6. 

m -m o 

= F(y -X ) 
s m 



V 



4 
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where u vector input to both system and model 

- vector output of the system 

- vector output of the model 

fg - vector structure of the system 

f - vector structure of the model 

— m 

V - scalar function to be minimized. 

A simplified version of the above v/ill be 



yg(t) = f(u(t^,t), t, £g) 



( 6 . 2 ) 



where the values of n are unknovm and the values of o 

— s 

will be determined by minimizing V, which is often of 
the form 

= _1 lCZs('C')-yj^(X)l'^d't (6.3) 

^"■^o '*^0 



and for optimum £ (best approximation) _6V = 0. The 

vector 2 thus found will only be optimum in the interval 

t^ < tr <t and m.ay not apply outside that interval. Also 

it may not be true for non-linear systems that 

lim J 2 = £ (6.4) 

"t- — ►oo rii s 

as one might expect from linear systems. The optim.ura 

value of 2 in Equation (6.3) could be found by a maximum 
m 

gradient seeking algorithm (first order gradient method) 
or by a second order gradient method ([29] 
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The second approach is to use linear estimation 
techniques, such as a Kalman filter to the linearized 
form of a modified state representation of the system. 

Considering a linear system driven by noise, Bryson 
and Ho [29, Chp. 11 and 12] present the following expressions 
for the estimate of the system: 

l) Linear Multistage Systems 



X = F.x. + G.w 
“i+1 “i“i 



X - state vector 



E[x 1 = X 
~o “o 

-X 1 [x -X 1 
o o 0 0 



= M 



( 6 . 5 ) 



W 



- environment noise ELw] = w 



E[v/ . “W . ] [w . -w . ] ^ 
1 - 1 ^ 0 “3 

_ _ T 

ePw.-w.I fx -X 1 
1 — 1 -* — o — o 






= 0 



F- G. - transition matrices 
—1 —1 

and z, = H. x, + v, 

-1 -1 “i “1 

z - measurem;ent vector 



- mieasurement matrix 
V, - measurement noise E[ v . ] = 0 



( 6 . 6 ) 



ECVi^l = 



E[x "X 1 v^ = 
o o 



0 



t 

i 



f 
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Assuming the sequences i of the different vectors 

to represent Gauss-Markov random sequences, they present 

the solution for the estimate of the state x x asi 

X. =x. +K.(z. -H;X.) 

—1 —1 — 1—1 — 1 — !■' 



Si+i = «i 



T -1 

K. = P.H. R. ^ 
— i —111 



( 6 . 7 ) 



Zi = + h.^r/^h.] 



-1 



1-1 



Mi - MiHihH.H^Hi'' + Rp-1 
T T 

K. , = F.P.F- + G.Q.G. 

— 1+1 — 1 — 1—1 — 1 — 1—1 

where P. is the error covariance matrix of x. after the 
—1 —1 

measurement while is the error covariance matrix 
before the measurement. Similarly x^ is the estimate of 
the system after the measurement while x^ is the estimate 
of the system before measurement. 

2 ) Linear Continuous Systems with Continuous 
Measurements 

The mathematical expression for this type of 
system is obtained by letting the time betv/een stages go 
to zero. Therefore, for the system 



X = Fx + Gw 
z = Hx + V 



( 6 . 8 ) 
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with F = lim F(t- +a » "t- ) - I 
A-0 7 — ^ 



G = lim G. 
A-*0 — 



( 6 . 8 ) 



Under the assumption that x represents a Gauss-Markov 
randon process, the expression for the Kalman filter 
will be 

^ + K(z - Hx ) 

T 1 

K = P H R 

p = fp + pf'I' + gqg'^-ph'^ S’^ h £ 

P <"=0) = Po 

If, for the system as specified in (1.11) to (1.13) 

- L (xiil»p') + iL* (1,11) 

£ = 0 (1.12) 

z = H X + V (1.13) 

the above technique is to be used, it is necessary by 
first to form the extended state vector 

X = 



r *1 
X 



(6.10) 



R 

The new state equations are as follows 
X = f(x»u) + w s[x(tQ)l = Xq 

^ = H x + V 

where in addition to 6.10 



( 6 . 11 ) 



w = 



w ' 
0 



( 6 . 12 ) 




V 



- 68 - 



f = ff*' 

.0 J ( 6 . 12 ) 

H = [H* o] 

Now the state equations are linearized in order to 

use Equation (6.9) or: 

• __ 

X = f(x»u) + w + K(^-H x);^("to) = 2Lq 

(6.13) 



K = P H^R"^ 



P= 6fP+P6f*^+Q-P H^R"^HP; P(t^) = P 
6x 6x 

The next simplification will he, considering that in 
practice most measurements are discrete in tim.e, to 
approximate the above differential equations by a set of 
difference equations and use the form.ulation for a linear 
multistage processes or, in other words, for the system 
below : 

= 2ii +A.f(2Li»lii) + 



(6.14) 



H x^ + V 



The resulting filter will be, assuming for convenience 

that 

E[A*Wjj^l = 0 

^i+1 = ^i ^*%Ui - ii ^0 

Hi = Eili^(H -^)**^ 



(6.15) 



4 
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P. = M. - A-M.h'^(H M. H?a+ R- a H M. 

1 J- 1 

M. - = 6 (x.+A'f(x. ,u. )) P. 6 (x.+A’f(x. ,u. )) Q*A 

“^+1 6 i “ 111 - -1 -1 



(6.15) 



where R and Q are constant matrices. If the value of A. 
is sufficiently small, this formulation v/ill represent 
closely Equations (6.13) ^.nd will be much easier compu- 
tationaly . 

In the above formulations, the partial derivatives 
m.ay be evaluated in a nominal path or if greater accuracy 

A 

is desired they are evaluated at x = 2£* 

Bryson and Ho [29l mention several authors v/ho have 
achieved success using Equations (6.13)and (6.15) for 
non-linear systems but at the sam.e time m.ake the warning 
that convergence m.ay not be obtained if the initial guess 
is poor or if the disturbances are so large that the 
linearization is inadequate to describe the system.. 

From the foregoing Equations (6,2) and (6.3), 
Equations (6.13) and Equations (6.15) i't v/as decided to 
use Equations (6.15) due mainly to computational 
convenience. Due also to the poorness of the available 
data, as mentioned in Chapter 5? it was decided to abandon 
the identification in heave and coupled heave and roll 



( 
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and concentrate exclusively on the identification of roll 
motion parameters whose data at a frequency of excitation 
of .66 cycles/second (roll natural undamped frequency) 
seemed reasonably good. 

As mentioned in Chapter 4, Equation (4,17) the state 
equations for roll motion, ignoring for simplicity the 
dependency of the system on the derivative of the input 
are: 



p = Kp p + [ - W.BG ] (^-fy)+ 



■A+I 


-K*+I 




P XX 


P XX 


p XX 


W . BG 
•p-*-Ixx 




1 p1p| + [1/3JK + 


-K* + I 


-K*+I ' 




P XX 


P XX 



(4.17) 





3 




P 



Non-dimensionalizing the above equations v/ith the 



frequency of excitation the following is obtained; 



P 




* l' 




■D 




K 


• 1 - 


W . BG 


1 






_-Kp+Ixx 




iO 




-K'+I 
P XX 


^ ■ 


-K* + I 

p XX 







W . EG 1 


<b 


l/2!Kp|p| 




-k'+i 
'^p -^xx 


1 





( 6 . 16 ) 
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1/3!K ^ 


1_ + 


1/31V/.BG 


1 


(‘J'-jI ■ 


1/31W.BG 1 


- 




(X>^ 


-K* + I 

p XX 




^-Kp+Ixx 



[iH^ 

Based on results of Appendices 1, 2 and 3, the 
following non-dimensional parameters will be, for 
= A . Iks sec . j 



K 

r 


1 




-.0042 


.01 


x(3) 


-V^xx 


CO 










K 

P 


1 




-.989 




x(4) 


-^^p-^^xx 












W . BG 


1 




2 . 744 


10 


x(5) 


-K' + I 

P XX 


u)2 










1/2 p 




= 


-.126 




x(6) 


-''p'^lxx 












1/3 IK ^ 


1 


= 


3.123 


10 


x(7) 


-S-^^xx 


— S’ 
CO 










Also and 


to 


have the values of 


all the state variables 


between ±1, 


it 


was chosen 


that 






* 

2_ 




.1 y(i) 








R 

cO 


=: 


.1 x(l) 




(6.18) 
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.1 y(2) 
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=.lx(2) 

(6.18) 

=-lu(l) 

The complete set will then have the form ; 

.ly(l) = .01x(3).lx(l) f [x(4)-10x(5)] [.lx(2)-.lu(l)] 4 

+ 10x(5).lx(2) + ) 

+ x(6).lx(l).l|x(l)| - [l0x(7)-.167 10 x( 55 
[.lx(2)-.lu(l)] .167 10x(5) [ilx(2)] ^ 

.ly(2) = .lx(l) 
or, rearranging 

y(l) = .01 x(3) x(l) + |x(4) - 10 x(5)] |^x(2) - u(l)| + 

+ 10 x(5) x(2) + (6.19) 

+ .1 x(6) x(l) ABS x(l)t[,l x(7) - .0167 x(5)] 
[x(2) - u(l)]^ - .0167 x(5) x(2)^ 

y(2) = x(l) 

Before using the data, this model v/as tested by si- 
mulating noisy measurements — only x(l) and x(2) were 
measurable in the simulation and the control u(l), 

with a format similar to the one obtained in the data 
collected. 

Also simplified forms of the state equations were 
tested. Theii* expression v/are : 

Non-linear equa t ions of mot i on to Identify linear terms 

y(l) = .01 x(3) x(l) + [x(4) 10 x(5)] |’x(2) - u(l)j + 

+ lo x(5) x(2) + .0126 x(l) ABS( x(l) ) + 

( 6 . 20 ) 

+ [.03123 - .0167 x(5)][x(2) - u(l)] * 

4 .0167 x(5) x(2)^ 
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y(2) = x(l) (6.20) 

Linear Equations of Motion 

y(l) = .01 x(3) x(l) + [x(4) -10 x(5)] [x(2)- 

u(D] + 10 x(5) x(2) (6.21) 

y(2) = x(l) 

Equations (6.19) » (6.20) and (6.21) were tested 
varying the following parameters: 

1. Initial State Estimate - The "true values" of 

the parameters (x(3) through x(7))were set in the simulation 
at the values given by Equation (6.1?). The initial 
estimate was set either zero or Equation (6.17). 

2. Initial State Estimate Error Covariance - This 
is a diagonal matrix of equal elements which v/ere set to 

0 , 2. and 5 . 

3. Order of Magnitude of the Elements of R Relative 
to the Elements of Q - Both the$ matrices are diagonal and 
of eo^ual elements. The ratio of the elements of R to the 
elements of Q was set at 1, 2, 5 and 10. 

4. Noise Level in Environmental Noise - w and in 
Measurement Noise - (r , Set- at 0 and 1 

Then the collected data was used in the identification 
procedure introducing small variations in the param.eters 
mentioned in 1, 2 and 3 above around. a reasonable set 
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obtained in the simulation to obtain the best possible 
convergence and also to test the effect of initial 
conditions - or in other v/ords, the necessary accuracy 
in the theoretical prediction of the values of the para- 
meters, to obtain convergence. Some of the results from 
the identification are presented in Appendix 5 while the 
discussion of the results is presented in Chapter 7. 

Summarizing in this chapter some of the available 
techniques for parameter identification were discussed. 
Attention was focused in those that permit application 
to non-linear dimamic system.s. From those, one was chosen 
that having perhaps less accuracy in the results, is easy 
to implement. The equations of motion for roll went 
then through the necessary modifications to be in agreement 
v/ith the format of such formulation. The structure of 
the model and the formulation v;ere tested by simulating 
the excitation and the measurements, and varying selected 
parameters to obtain a better insight of what to expect 
from the identification. 
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CHAPTER 7 



DISCUSSION OF RESULTS 
1. Simulation 

The first equation to be tested v/as Equation 
6.19* The noise level was set to zero, the ratio of ele- 
KBnts of R to elements of Q, from now on expressed formally 
as r/q by similarity with the scalar case, equal to 1, 
the initial state estimate error covariance matrix at 
[ 2 ] and the initial state estimate Xq at zero (the ’’true 
value" was set in all simulations at the values expressed 
in equation (6*17) )♦ The values of the ostimabe blow up 
and after 200 measurements the computer indicated overflow. 
This type of behaviour with no noise decided the use of 
Equations (6.20) or (6*21). These two models wer^ both 
tested at no noise. Equations (6.20) wore first tested 
with x(0) " Eq.(6.17), r/q = 1 and P q = [^] « error on 

the estimate ( XHAT - XX ) where XHAT is the estimate and 
XX is the "true value" of the state , remained at zero during 
the 7/hole experiment. However P(3,3) (variance of x(3) - 
linear damping coefficient) kept increasing v/ith small 
oscillations. The other elements of P decreased quickly 
to less than .1 . The same behaviour was observed for 
Equations (6.21), showing that the influence of the 
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non-linear terms is small, at least for the no noise 
condition. 

The last experiment in the no noise series was con- 
ducted \7ith Equations (6.21) Pq = [2] , r/q - 1 and 
( XHAT )^= [01. 

The error ( XHAT - XX ) for all variables converged 
rapidly to zero - in less than 100 measurements V7ith the 
exception of (XHAT(3) - XX(3)) v/hose value remained 
constant throughout. Again the term P(3,3) steadily 
increased; v/hile all the other variances came to values, 
after 500 measurements, less than .060 (initial value •was 
2,0) the final value of P(3,3) v/as 2,202. The foregoing 
suggests that when no noise is present the system converges 
independently of initial conditions, v/ith the exception 
of Equations (6.19) and v/ith the exception of the linear 
damping coefficient x(3) v/hose identification seems very 
difficult to obtain. 

The rest of •the experiments v/ere performed v/ith noise 
level equal to 1 which is believed to be quite a high 
level due to the scattering obsev'ved in the measurement 
vector z. 

The first simulations v/ere taken v/ith r/q = 1 
XHATq = [0] and [3] and = [5J . Both simulations blew 



up at about the 140th measurenBnt. At that time p(3,3). 
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P(4,4) and P(5,5) were quite large, but different from 
one case to the other. The main point to be observed here 
is the conjugate effect of noise and not accurate initial 
e stimate . 

The next experiment put together Equations (6.20) and 
(6.21). For both systems r/q=r 1, [o] and ( XHAT )^ = 

(XX) ^ » Since initially the error covariance matrix v/as 
set to zero, it increased and then oscillated about 
values equivalent to the ones found in the no noise 
condition (less than .050), with the exception of P(3,3) 
which steadily increased to a final value at the 500th 
measurement of .249 in both cases, with tendency to 
incr’ease. The errors on the variables oscillates around 
zero with the exception of (XJiAT(3) - XX(3)_) v/hich remai- 
ned constant at about .035. Again it seems that the 
identification of x(3) v/ill be dificult to obtain. 

Next an experiment was performed v/ith Eq. (6.20) 
with r/q= 1, ( XHAT )q = i^) ^ but P^- [s] which diver- 

ged again confirming the influence of the initial con- • 
dltions when combined with noise. 

The next set of experiments wore performed with 
the noise level at 1 , as before, but r/q =: 2, 5 and 
10 . 

Comparing the results from Equations (6.2l) for' 
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r/q 1 and r/q 2, ( XHAT )^ _ [o] and in the later 

case the estimate does not blow up although the error 
on the estimate and the variance both increase for x(3) 
and x(4), v/hile for x(5) the error on the estimate and 
the variance decreases, v/hich demonstrates some impro- 
vement. The cases r/q = 5 r/q = 10 represent no advan- 
tage over r/q= 2, presenting the same trends. 

Similar simulation was performed v^ith Equations 
(6.20). Abain an improvement is noted when r/q goes 
up to tv/o but no further improvement is obtained v<hen 
r/q = 5 or 10. On the contrary, the Mecre asing in the' ' 
error seems to take place more slowly in the later 
case. In the case r/q = 2 the advantage of considering 
the non-linear terns when identifying the lineai^ para- 
meters is noted because P(3,3), which increased in the 
case of Equations (6.21), is seen to oscilate between 
approximately 2.02 and 2.16, and P(4,4), after increa- 
sing for the first couple of measm’emsnt, markedly 
decreases. The foregoing results were then compared 
with the case of P^= [5] , v/ith clear advantage for the 
former . 

The foregoing seems to indicate that the values 
r/q = 2 and ai’C the befit suited among all 

other combinations tested. At this stage small varia- 
tions could be tried to further improve the convei’gonco . 

But since the carac teris tics of the data are not 
expected to follow very closely the carac teristics of 
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the simulation, it seems a better idea to try to iden- 
tify the linear parameters with the data obtained from 
the experiments and then vary the values of R, £ and ^ 
in order to obtain a better convergence. But before 
that, one last simulation was tried for the same values 
of Q, R and P^, but for values of ( XHAT ) ^ closer to the 
true values : (XHAT) = fO, 0, -.2, -.5, .151 ( about 

Q U 

1/2 of the true values ). The number of measurements 
v/as also increased to 1000. It was noted that : 

i) P(5,5) decreased rapidly to a value of .4 
(after 10 measurements ) and then decreased oscillating 

until it reached a value of .057, The error on the esti- 
mate, after 1000 measurements, was oscillating between 
± .1 . 

ii) P(4,4) decreased less rapidly to a value of 
.566 ( after 30 measurements ) and then jumped to a 

value of .06 v/here it kept oscillating for about 350 
me asureraents decreasing finally to a value of ,025 where 
it stayed oscillating until the end. The error on the 
estimate oscillated between ± 1 for 400 measurements 
and then decreased to zero oscillating about that 
point , 

ill) P(3,3) increased to a value of 2.22 after 
400 measurements, oscillated for a while and finally 
decreased to reach at the end a value of 1.29, with 
tendency for further decrease. The error on the estimate 




i'**! 




f 




f 



■I 






I 

I 




-80- 



Increased to a value of 3.3 after 420 measurements* 
oscillated about that point for a while and finally 
decreased, oscillating to a value of 1 at the end. 

Note that the turnover point seems to be the 
400th measurement and that a set of initial conditions 
on the estimate closer to the true value really impro- 
ves convergence. This fact confirms the believe that, 
contrary to the linear case, the initial conditions on 
both the estimate and its covariance play an Important 
in the success or failure of the Identification. 

2 . Identification Using Experimental Data 

Follovdng the conclusions of last section 
a run was made with Q = , S— * 

The set of initial conditions was as prescribed by 
Equations (6.17). Simultaneously, runs v/ith the values 
of Q, R and P as Indicated below were also made. The 

— — — ’O 

set of Initial conditions, however, v/as not changed at 
this stage. It seemed preferible, since the above set 
is believed to bo, among the available predictions, the 
closest to the true value, to find first the best sot 
of Q, R, and changing (XHAT)^ afterwards. 

Values of Q, R and Pq : 

2 = [-ij 

R = [.1] , [.2] , [.3] 

[1-] , [2-J , [5-] 

in the following combinations of R and ; 
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( [.05], [2.] ) 

( ^ ^ ^ ^ ’ C® *J ^ 

( C.3 ],[!.] ) 

For the group ( f.2 ] ,[s,] ) the results obtained 

with the first measurements v/ere more or less as 
expected* P(5,5) decreased rapidly, P(4,4) decreased 
less rapidly and P(3,3) increased steadily* Then P(4,4) 
came close to zero and be c am e ne ga t ive for a while, then 
crossed zero again and went up to 2*48, decreased to 
*270 and increased again to the final value of *421, 

Also at the end P(5,5) = .059 and P(3,3) = 2.473 * 

On the other runs a similar behaviour was noted. For 
R- [. 2 ] and ?.= r^-] , P(4,4) crossed zero and kept 
decreasing up to a final value of -1.589* P(5,5) 
decreased to reach the value of .064 ( at the same 
point at v/hich P(4,4) became negative ) and then 
increased to reach a value of .123 at the end* P(3,3) 
increased steadily to reach 1.462 at the end. For the 
case R = [.2] , P^=r [5,] it was oven worse since at 
the end P(3,3) = -43,225 and XHAT(3) = 252 ,70 .* .For 
R;ir|.l[ , P^=|5,| , P(3,3) increased steadily, P(4,4) 

decreased to a valvio of ,173 and then increased again 
to reach at the end .425 with XHAT =-14.513 ( ? ). 

P(5,5) decreased to a value of ,028 and kept oscilla - 
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ting about that value, but with XHAT(5) alv;ays Increa- 
sing. For 5.= behaviour was similar 

but not so extreme due to the smaller value of P • 

— o 

For R z: * P(3,3) always increased 

while P(4,4) and P(5,5) decreased first to low values 
( .60 and .20 ) and increased afterwards, never becoming 
negative. For the last run of this group { Rt=j^.3^ , 

)» the type of behaviour was similar to 
( R=f.2] , ^=Ql.]] ) but much more extreme - — the 

last recorded value of P(4,4) v/as -226,1 at the 350th 
measurement because some cycles later the computer 
stopped by overflow. 

The fact that P(4,4) became negative posed the 
question of the correctness of the computerized scheme. 
HoT/ever the program after thourough checking have been 
v/orking successfully with linear systems. Besides no 
such strange behaviour Y/as found during the simulation. 

Therefore the program v/as checked again in detail 
but nothing could be found able to explain the fore?-".n‘. 
going. 

It v/as sugested then that as P(4,4) approaches 
zero, the computer could introduce some errors thus 
making it negative. The process would bo self sustaining 
afterwards. Assuming that to be the cause, the solution 
proposed v/as, v/henevor P(I,I) became negative, to 
substitute it by its positive counterpart. Accordingly 
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the following statements were introduced in the Kalman 
filter ( KLMN subroutine ) : 

DO 21 I 1, NX 

IF( P(I,I),LT.0.0) P(I,I)= -P(I,I) 

21 CONTINUE 

The effects of this procedure were then verified 
on the follov/ing combinations of R and ^ : 

( [-2] . [l-J ) 

^ C*®! * ^ 

( [.3] . [' X.], ) 

For the second case, this had as a consequence to 

hold XHAT(4) within bbunds longer than before, only to 
diverge more violently towards the end ( XHAT(4) == 735« ). 
The other tv/o parameters had also large values* For the 
first case, the influence of the above procedure is more 
evident, since originally XHAT(4) decreased, crossed 
zero at approximately the 100th measurement ( point 
where P(4,4) became negative ) and v/ent on decreasing, 
With the presence of the above statements, it behaved 
in the same v;ay for the first 100 measurements, then 
stayed at a value of .8 ’■’■ntil the 300th measurement, 
starting a slow decrease through zero to reach the 
value of -2.605 at the end . Acordingly the variance 
P(4,4), after touching zero increased to a final value 
of .242 . The sarnie type of behaviour is also noted in 

I 
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both XHAT(3) — diverging v/ith negative values before, 
with P(3,3) increasing and being constant now with 
P(3,3) increasing also, and XHAT(5). - — going av;ay 
from zero with positive values before and going away 
from zero but v/ith negative values and at a much slower 
rate, with the introduction of the above statements, 
P(4,4) seemed to be bouncing between zero and some 
finite positive values, but since the values of P were 
only recorded every 50 measurements, no valid conclusion 
can be taken about that fact. P(3,3) alv/ays increased 
and P(5,5) remained at approximately .067 • 

It v/as decided then to modify the structut'o of 
the state equations to check if another form of the 
equations v/ould converge to a solution. In Equation 
(6.16), for example the tvjo groups forming the 
coefficient of x(2) - u(l) could be considered as 
one parameter, being the new fox’ms of x(4) and x(5) ; 

• • K - Vi.EG 1 

10 X ( 4 ) — 

"^p + ^xx 

W'. BG 

10 x(5) = 

^xx 

And the new form of the state equations would then be : 
y(l) = .01 x(3) x(l) -f x(4) |^x(2) - u(l)J 

+ 10 x(5) x(2) + .0126 x(l) AES( x(l) ) 

^ 3 (7.2) 

+ [ .0312 -f ,0167 x(3)] Qx( 2) - u(l)] 

- .0167 x(5) x(2)^ 



to 



1 






(7.1) 
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y(2) x(l) 

For comparison only the following valued of R and 

were considered ; 

— o 

1) R = and ^ [5.] 

li) R=[.3j and I'l.'] 

Also and to check the behaviour of the system at 
high values of Q and R, the following run was made : 

ili) Q=[1.J , R=[2.J , Pq.z[1,] 

The main caracteris tics of these runs was that the 
parameters did not oscillate, or oscillated v/ith a very 
long period. Only in the case R=Q.2j , , 

XHAT(3) v;ent over 1. at the 310th measurement, 
approximately. XHAT(4) and XHAT(5) which started at 
-.373 and +.274 respectively, crossed each other and 
zero at the 520th measurement reaching at the end the 
values +.359 and -.422 and showing in between a 
behaviour completely linear. P(3,3) increased to a value 
of 5.034 at the 100th measurement and then decreased as 
steadily as it had increased to a final value of 3.412. 
Similarly, P(4,4) increased to a value of 5.003 in the 
first 15 measurements and then decreased to a final 
value of 4.838. P(5,5) from the st^art decreased oscilla- 
ting until it reached a final value of ,304 « This 
looked like being the object of the search apart from 
the fact that the values the parameters reached at the 
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ond were not to be expected on physical grounds, 
specially XHAT(5) ( this 7/ould mean that the center of 
bouyancy had to be above the center of gravity which is 
not true from the calculated curves ). A. possibility to 
consider is that the variables go through a cycle of 
very long period and therefore the convergence v/ill take 
place later in the process. For the case R , 

^ [l.J , only the value of P(5,5) decreased and the 
parameters changed their value very slowly and therefore 
nothing could be said about such a process. It could be 
an oscillatory motion of a period even longer than 
before or it could be a steady, but very slov/ly divergent 
process. In the last case considered, the pai’ameters did 
not change but the diagonal elements of the covariance 
matrix' Increased and quite rapidly to reach values of 
[1.91 , 2.239 , 48.289 , 48.195 , 32.620] . Mote that 

that the variance of XHAT(l) and XHAT(2) was, for all 
the cases considered, simulation inclixded, of the order 
of .025 and .012 after the first 20 to 30 measurements. 
Therefore this last case is of no consequence <> 

In sumriiary in this chap'ter the resultvs from the 
identification were presented, both for the simulation 
carried out and for use of the data. It y/as pointed out 
of the difficulty of the system to converge. A set of 
values of the different covariance matrices was selected 
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from the simulation studies as being those that offered 
more possibilities iVhen using real data. Hovrever the 
results didn’t meet the expectations, being necessary 
to proceed in looking a set that v/ould make the system 
converge • In stead it v/as found that som.e of the 
diagonal terms of ^ would become negative. First it was 
thought to be resultant from a mistake in the compute- 
rised scheme and later to errors Introduced by the 
computer itself and corrected accordingly. Even then 
no solution could be found after 1000 measurements . 

Then another form of the state equations was obtained 
by combining two of the parameters. This later system of 
equations presented som^ interesting carac teristics, 
but nothing could be said in definite. 



I 
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CHAPTER 8 

CONCLUSIONS AND RECOfilEITOATIONS 

From Chapter 7 it was quite clear that convergence 
of the estimate of the system to its true value was, in 
general, not possible to obtain. Reasons that could 
account for this failure in convergence were then pointed 
out as follows : 

1 - Possible mistakes in the computerized 

scheme . 

2 - Not accurate prediction of the initial 
Condi tionsc of vthe ' system. 

3 - Data that measured the dynamical behaviour 
v/as not accui'*ate enough. 

4 - The input to the model didn’t really haVe 
the same caracteristics as the input to the system 

in the experiments. 

To those reasons above anothei’ may be pointed 
out that could be partly responsible for the non- 
convergence of the system, namely : 

5 - Incompatibility between the system 
structure and the particular technique that was chosen . 

V/lth respect to' the first reason, the fact , 
mentioned already that this same procedure was used 
before in the identification of linear static systems 
successfully, although not a suficient reason to 
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completely exclude this possibility, at least decreases 
the probability of such an event. Therefore it is recom- 
mended that if future use is intended to such a formu- 
lation, the program be rechecked for mistakes. 

The not accurate prediction of the initial condi- 
tions can be a reason of some standing if the methods 
used for prediction fall under the category of empirical 
or semiempirical formulations from which the final 
result have not the reliability required. For the calcu- 
lation of the motion parameters for roll, for example, 
the three categories in which Chapter 3 was divided 
also provide a good basis for classification of 
accuracy, starting with the hydrostatic coefficients 
at the top and the non-linear coefficients at the 
bottom of the scale. It is believed that for the case 
of excitation due to regular waves a formulation as 
presented by Ref. IS and 19 for linear added mass, 
added moment of inertia and damping coefficient is 
completely adequate in the sense that the system will 
converge if the values of those parameters constitute 
the only approximation involved. This c inclusion is 
reached with Information on the simulations carried 
out in Chapter 7. 

V.hien trying to digitize the data from the towing 



tank experiments, it v/as believed that this item 
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played a major role in the success or failure of the 
identification. This is not so anymore, when one compares 
the data from the tank with the measurements generated 
in the simulation, where the ” level "‘of noise is much 
higher than in the former case. This is not the same as 
saying that it is not necessary to care about the 
experiments. If nothing else a careful planning will 
save a lot of time and labor, waisted otherwise trying 
to correct past mistakes. 

The accurate determination of the input to the 
system as a function of time seems to play an important 
role in the success of the identification. Intuitively 
if the behaviour of tv/o systems' ist to be compared, it 
is necessary that the input be the same. Unfortunately 
this was not the case. Both methods presented in 
Chapters 3 and 4 ( using the Haskind Relations or 
assuming the major part of the excitation to be due to 
changes in bouyancy :due 1 to : the ..passage of the v/ave ), 
will not give very accurate results, the first because 
it does not present the instaneous values of the input, 
but an expression that assumes before hand an harmonic 
tims dependency, the second because, although it 
expresses the input as the instantaneous function of 
a quaiitity being jneasured in the experiment, exprimes 
a very rough approximation of such dependence. However 
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the second was used in the identification procedure and 
it is still believed that, por this particular case the 
second yelds better results. 

The incompatibility between the system structure 
and the method of identification exists because some of 
the assumptions set forth in the derivation of this 
particular method of identification were not met by the 
system. One of such assumptions was that the noise in 
the system, both environmental and measurement noise 
wore uncorrelated with each other and with the state 
of the system ( Equations 6,5 and 6.6 ), This last 
point v/as not met by the vehicle structure since, for 
example, finite water depth effects, interference 
betv/een the vehicle and the v/alls, etc., would in one 
v;ay or another be correlated v/ith motions of the 
vehicle. This fact as a reason for failure to achieve 
cdnvergence gained Importance when comparing the 
results of the slrmalation with the results of the 
identification, where the former, being bad, were a 
bit more encouraging than the last. 

With the above in mind, it seems that in order to 
improve the application of the identification techniques 
to problems of the type discussed herein, it v/ill be 
necessary : 

1 - Generate a vehicle environment that can be 
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better approximated as an uncorrelated random process* 
This sugests the possibility of using random seas to 
generate the data for the identification. 

2 - Use another method of system identification, 
such as the model ref fere nco technique briefly 
described in Chapter 6, which do not make such restrict 
tive assumptions. 
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APPEKDIX 1 

HYDROSTATIC COEFFICIENTS 

1 - DISCRIPTION OF THE COMPUTER PROGRAMS TO OB- 
TAIN A POLYNOMIAL SURFACE REPRESENTATION 

2 - VALUES OF THE HYDROSTATIC COEFFICIENTS FOR A 
MODEL OF THE ” MARIl'JER " HULL FORM 
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Figure Al.l 

Hydrostatic force vs, angle of rotation and vertical position 
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VALUES OF THE HYDROSTATIC COEFFICIENTS 
Since the identification of the motion parameters in 
coupled heave and roll v/as not performed, it was necessary 
to calcula.^ the hydrostatic coupling coefficients. To cal- 
culate only the coefficients Z^, , and K^3 a much 

simpler method than resorting to the Kerv/in polynomial re- 
presentation and that will probabilly yeld more accurate 
results can be used • The reason for such is that in com- 
puterized solutions using least squares, it is not usual- 
ly possible to specify some a priori knov/ledge that one 
might have about the surface or the curve. For example in 



roll K(iz5=0) and K^2 



symmetry. Also 






are zero due to port and starboard' 



** . A least square method gave 

nonzero values for the first tv;o coeff iciente s. 

Looking just for the secon non zero term in the Tay- 
lor's expansion and knov/ing the first it’s much easier to 
calculate manually the missing term. Using this procedure 
and being, for roll - W-x-GI.' -2.159 LB-:c-FT, the value that 

will approximate the mioment curve better will be ; 

Between ^^1.175 l/3'-::-Ko3 ^ 6,814 LB-:;-FT 

^^1,350 ^ -1.254 LB-”-FT 

^z-t,525 ^ -3.007 LB;:-FT 

Since the roll angle will go over ,175 but does not 
reach .350 in the Gxoerimonts the value choose n was: 



1/3 'Ko3 - 6.814 LB-::-FT 

Similarly for heave the values found were 
Zz = - f = -194.709 LE/FT 

1/2' 3^2 -84.528 LB/FT^ 
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APPENDIX 2 

LINEAR HYDRODYNAMIC COEFFICIENTS 

1 - DISCRIPTION OF THE COMPUTER PROGRAMS TO CAL- 
CULATE ADDED MASS, ADDED MOMENT OF INERTIA AND DAI.^PING 
COEFFICIENTS FOR T\'iO DIMENSIONAL CYLINDERS IN HEAVE, SWAY 
AND ROLL 

2 - VALUES OF ADDED MASS, ADDED MOMENT OF INERTIA 
AND DAl^PING COEFFICIENTS FOR A MODEL OF THE ” MARINER " 



HULL FORM IN HEAVE AND ROLL 



t 



1 
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PROGRAJ.l INPUT 

Cards 1-4 FORMAT (20A4) 

The purpose of these cards is to provide titles for 
the computed output. V.lien titles are not desired blank 
cards should be provided. 

Card 5 -NOK,MAXB,NH F0RHAT(3I6) 

NOK - Total number of frequencies, at v/hich the cal- 
culations will be performed. (MAX 50) 

MAXB -Integer correspondent to the horizontal plane 
of maximum beam. Start the count at the keel, 

NH - Number of depths at which the calculations will 
be performed, (MAX 10) 

Card 6 - CAY(K),K 1,N0K FORMAT (5F12. 7) 

CAY - Nondimensional frequencies - a/g 

- radian frequency 

a - half the beam for a floating cylinder 

- half the draft for a submerged cylinder 
g - gravitational acceleration 

Card 7 -(DSP(J),J 1,NH),CR F0RMAT(5F12.7) 

DEP - Values of the depths of submergence, measured 
fro^i the svirfaco to the top of the cylinder. For a floa- 
ting cylinder, NH 1 and DEP 0 , 

Card 8 - MD,K0N FORMAT (21 6) 

MD - Control integer - 1 for a submerged body 

2 for a floating body 
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NON - Number of segments in v;hich the right hand side 
half contour is divided (MAX 45) 

CR - Positive vertical distance from the center of 
rotation to the free surface, for a floating cylinder or 
from the center of rotation tb the intersection of the 
section contour with the plane of symmetry( on the side 
closer to the surface* 

Card 9 - X(J),J 1, (NON 1) F0RMAT(5F12.7) 

X(J) - Horizontal offsets of the end points of the 
line segments, measured from the vertical planr of symme- 
try to the right hand side of the contour, starting at the 
keel. The last point is the intersection of the section 
contour either v/ith the free surface ( floating cylinder ) 
or with the vertical plans? of symmetry( submerged cylinder ) 
Card 10 - Y(J),J 1, (NON 1) F0RMAT(5F12.7) 

Y(J) - Negative vertical offsets of the end points of 
the line segments, measured from th free surface, for a 
floating cylinder or from the intersection of the contour 
with the vertical plane of symmetry, for a submerged cy- 
linder, starting from the keel* Again only the right hand 
aide of the contour v/ill be considered* 

Card 11 - MODE F0RMAT(I6) 

MODE - Control integer - 1 - heave 

2 - sway 

3 - roll 

Card 12 - FORMAT ( 20 A4) 



I 
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Provides a title for the output of the program. If 
the title is not desired, a blank card should be provided. 

If more than one mode is desired, adicional pairs of 
cards 11 and 12 should be provided, accordingly. 

Card 13 - MODS FORMAT (16) 

MODE - Control variable 
0 - End of data 

Negative integer - A new geometry follov/s 
starting Card 5. 



I 







THL l“uLLOWii\iG PRuGRAiMS WcKt; IAKt:K rKUfs 
USiluC RuPJkT 3‘j51 

‘‘ NUMcKICAl lALCULATIGN UF AJiJEu MASS AKU DAMPING 
CCJc.FFIC IlNTS uF CYLIKDcKo JSCILLATIKG Ii\ or DElOW 
A Fh L c SuRpACl ** 
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APPEITOIX 3 



NON-LUfEAR MOTION PARAI4ETERS 



1 - DETERMINATION OF HYDRODYNAMIC COEFFICIENTS 



USING EMPIRICAL IffiTHODS 



-146 



1 - ADDED MSS AND ADDED MOMENT OF INERTIA 

Using Ref* 20 p* 61 to 63 - Heave 



U 

A m = --- y^(dx) 



(A3.1) 



where ^ - water density 1.94 LB'-'SEC^/PT"^ 

L - lenght of the model 5.5 FT 

- coefficient from fig. 17 ^ 1. 

y - contour of the waterline, aproximated by a 
parabola of the form 



y- B/2 



1 - 



X 

L/2 



n 



(A3. 2) 



,79 FT 



B -•* beam of the model 

n - C„/(l-0^) 

- Y/aterplane coefficient - Ay^j^/B-::-L = .724 

- coefficient, calculated according to 



= 1/ [ 1 +.8-B/L] = .895 (A3.3) 



the following form may bo obtained 

A m = o 772-:;-^,-;:- jO -;:-L--B^-::-C^^/( 1 ) ( A3 . 4 ) 



Plugging in the values 

Am = 1.41 LB::-SScVfT 

2 

Being the mass of tho vehicle - m 1.43 LB-;i-SEC /FT 
this x’osult may be chocked by measuring the natural period 
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of hoaving motions and using results from the hydrostatic 
calculations 2^ - ~ “ -194.05 LB/FT 

i 

8.29 SEC(-l) 

The value found for the natural undamped frequency v/as 
8.038 (SEC-1), 

Similarly for roll, the added moDisnt of Inertia is 
expressed as a percentage of the moment of inertia. This 
is found, from data ref faring to full scale vehicles, acor- 
ding to the formula 

^xx " V^/g-::-((B^-::-Cw)/ll.4-::-Cij'^ ) TON-K-SEG^ (A3. 5) 

= .02855 LBSEC^ 

From fig. 24 of the above refference 

A - . 286 -'>-I^j^ 

Therefore 

I + Al = .03665 LB-;:-SEC^ 

XX XX 

From the hydrostatic calculations 
K^=-Vr::-GM = -2.159 LB-::-FT 
The resulting natural frequency is then 
On "" 7.7 SEC(-l) 

while the value found in the measurements was 




m +A m 



lo w — 4.14 



SSC(-l) 




e 



I 




.1 
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The values of added mass and of addes moment of 
inertia, calculated in Appendix 2 were : 

A (3, 3) = 1.24 LB-::-SSC^/FT 

A (4, 4) = .0176 LB-::-SEC^i;-FT 

The value of roll added moment of inertia is seen to 
be quite different the results of both m.ethods are compared. 

2 - DAMPING COEPT'ICIENT 

For heave the following expression, taken from page 
163 of the above ref fere nee was used : 

-l/2j0 2'f7lv(L.B.C^)?/(TA);('t.LA) (A3. 6) 

where the parameters are as before and also : 

- vjave lenght = 3.11 FT for u;— 8.04 SEC(-l) 

- coefficient taken from (20) fig. 52 =: .460 

T - draft of the vehicle = .28 FT 

lV, 4' coefficient taken from (20)fig. 53 m .200 
For the values above : -2.336 LB-::-SEC/FT 

while v/hen calculated by the methods of Appendix 2 

Z^ = -7.4543 LB-::-SEC/FT which is three times the 
value presented above. 

For roll, the calculation of linear and non-linear 
damping coefficients were carried simultaneously, for each 
transverse section of the vehicle, using Equations 21, 22, 
23 or 24 of page 148 of the above refference, deppending on 
the shape of the section.. The value of the damping for the 
whole vehicle v/as obtained by summing the sectional results 



1 
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over the lenght of the vehicle in a " strip theory ” way* 
Therefore, using the following equations the damping is 
calculated. 



Kpi -l/2^b^Al^ 



, I , g(Bi/Ti-2.8)+4 
oi o, 1 ill (b /g)t 10 ^ 



(A3. 7) 



^Di 



Si 



72(B./T.-3.07f+-ll 



36(BVT.-2. )' 

Si < >^0 ’ Vh > I P K r~,-T5rS-r 



(bi/s) .‘cJ w 



(A3. 8) 



(A3. 9) 



-^1-1 



Si 



.57|p| + 



g(Bj^/T^-2.82) 4 



(bi7s)T 10^ 



(A3. 10) 



a 




where Eq. A3. 7 and A3. 8 are used for rectangular sections 
of different value of * Eq. A3. 9 is used with 

elliptical sections and Eq A3. 10 with terminal section of 
shape close to triangular, the following values were 
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found for linear and non linear damping coefficients: 



Kp -.00786 LB-::-FT-;:-SEC 

= -.06835 LB-;;-FT-;5-SEC 

= -.01605 LB-::-FT-;:-SEC^ 
P|P| 

The values found in Apenndix 
coefficients v/ere respectively -. 



- O) - 4.145. SEC(-l) 

- cu = 8.290. SEC(-l) 
all frequencies 

3 for the linear damping 
0022 LB-:c-FT---SEC' 

0363 LB-rcFT-::-SEC 



For the nonlinear heave damping coefficient, the va^ 
lues for each section v.'o:>e obtained by considering each se- 
ction substituted by an equivalent cylinder of the same 
diameter v/hen the model moves av;ay from the water ( assu- 
ming the only drag present to be separation drag ) and sub- 
stituted by a cylinder of the same lenght of wetted con- 
tour when the model moves inte the water ( assuming the on- 
ly drag present to be skin friction drag ) • Data from 
infinite cylinders was then used, the correspondent values 
of drag coefficient averaged over a cycle, and the results 
summed over the lenght of the model. The velocity entered 
in the calculation of the Reynolds Number was a RMS type 
of velocity taken from the records of the experlmexits. 

The final value obtained was: 

= -1.767 LB-::-SEcVfT^ 
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APPENDIX 4 

CIRCUITRY AND DIGITIZING PROGRAM 

1 - PROGRAM TO DIGITIZE DATA AND CORRESPONDENT 

PATCHING IN THE ANALOG COMPUTER 
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APPENDIX 5 

IDENTIFICATION PROCEDURE 

1 - DISCRIPTION OP THE COMPUTER PPOGRAI^S USED IN 
THE IDENTIFICATION PROCEDURE 

2 - EXAJvlPLE OF THE BEHAVIOUR OP VARIABLES AND RE- 



LATED VARIANCES DURING IDENTIFICATION 
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